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COMBINING OF PROFESSIONAL GROUPS 


HE Professional Group Task Force recently ap- 
BF vointen to survey the desirability and feasibility 

of combining Professional Groups have, after a 
study, come to the following findings: 

1. There is no significant general trend among Profes- 

sional Groups to combine. It appears that the scopes of 
most of the Professional Groups are sufficiently distinct. 
It is felt that it is highly desirable that all Professional 
Groups consider carefully the possibility of joint meet- 
ings wherever and whenever feasible. 
_ 2. There is no significant general trend among Profes- 
sional Groups to combine TRANSACTIONS. Several groups 
have discussed combining certain issues because of pos- 
sible interest to the joint membership. 

3. In order to assist and also perhaps to guide individ- 

ual Professional Groups in planning joint activities, it is 
recommended that Headquarters prepare a tabulation 
listing both horizontally and vertically all individual 
Professional Groups and placing at the intersection the 
numbers of IRE members belonging to both groups. This 
should convey quickly information about split of inter- 
est between two Professional Groups and might point in 
the direction of joint efforts on a quantitative basis. 
_ 4. It is recommended that each Professional Group 
include in its own TRANSACTIONS those abstracts and 
papers appearing in other TRANSACTIONS which may be 
of interest to its membership. 


‘ LIAISON WITH PROCEEDINGS 


_ As stated in the last issue of the TRANSACTIONS, Roy 
Spencer is serving as liaison between our quarterly and 
the PROCEEDINGS. When a technical paper in our field 
is submitted to the PROCEEDINGS, he is one of the three 
expert readers appointed. In this fashion, PGAP is rep- 
resented, we are able to keep informed of papers submit- 
ted, and we have a voice in deciding where they should 
be published. This arrangement has thus far been set 
up only for PGAP. 
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The Administrative Committee has been concerned 
for some time with the subscription rate for the TRANS- 
ACTIONS for other than PGAP members. Any arrange- 
ment that is made requires complete approval of IRE 
headquarters. It appears as if a very considerable num- 
ber of subscriptions can be obtained, and the rate estab- 
lished can determine a substantial amount of income. It 
is hoped that more information will be available by the 
time of the next issue. 

More local chapters are desirable. San Francisco, New 
York, Long Island, New Jersey, New England, and oth- 
er areas have been remiss in getting some activity start- 
ed. The Administrative Committee would welcome pe- 
titions for local chapters but each area seems to require 
some special spark to get under way. J. S. Brown of An- 
drew Corporation is responsible for local chapters; en- 
terprising organizers are invited to write him for further 
information. 

More Institutional Listings for the TRANSACTIONS are 
desirable. These are the company names and their 
spheres of interest as listed on the back cover. This is 
the only approved advertising we can employ, but at the 
present rate of $25.00 a listing per issue and $75.00 for 
the year, the possible return from this kind of advertis- 
ing is well worth going after. Our circulation certainly 
justifies the small expense on the part of the advertising 
company, particularly in the light of the well-defined 
audience exposed to the listings. Robert Helliwell of 
Stanford University is in charge of these listings and he 
should be contacted for information or arrangements. 

WESCON-1954 has made $800.00 available to PGAP 
for publication of papers presented at the Los Angeles 
convention this past summer. Papers will be subjected 
to our standard screening procedures and if accepted a 
WESCON acknowledgement will be given in a footnote. 
If we are too meticulous and publish too few WESCON 
papers the unspent balance of the $800.00 allotment 
must be returned. Even when the cause is just, it is al- 
ways sad to contemplate returning money. 


NATIONAL MEETINGS 


The Professional Group on Antennas and Propaga- 
tion will sponsor four sessions at the 1955 Spring Na- 
tional IRE Convention in New York. Two of the ses- 
sions will consist of five papers each on Antennas and a 
third session will consist of five Propagation papers. The 
fourth session will be a panel discussion of “Extended 
Range VHF and UHF Propagation,” chaired by Jerome 
Weisner. It is hoped that panel members will be re- 
cruited from Stanford, Bell Labs., the National Bureau 
of Standards, Cornell, and M.I.T. The availability to 
PGAP of four sessions at the National Convention 
should be very gratifying to our membership as a tribute 
to the high activity of our Professional Group. 

Dr. Bolljahn of the PGAP Administrative Committee 
is in charge of all meetings, symposia, etc. We have been 
advised that he has prepared a list of about ten national 
annual meetings that PGAP might conceivably sponsor 
either individually or jointly. It is hoped that his list 
will be available shortly. 

After the Spring IRE, the next important PGAP 
meeting will be the Spring URSI in Washington which 
we have recently been invited to once again co-sponsor. 


GROUP CHAPTERS AND CHAPTER NEWS 


The number of PGAP Chapters will soon be increased 
by 100 per cent with the addition of Chapters at Phila- 
delphia and Washington, D. C., and Los Angeles Orange 
Belt sub-section. 

The Philadelphia Chapter has received the approval 
of the PGAP Executive Committee and will start its 
program of activities in the near future. Those members 
in the Philadelphia area desiring more information on 
the Chapter activities should contact the organizer Mr. 
Charles Polk of the Moore School of Engineering of the 
University of Pennsylvania. 

The Washington, D. C., Chapter has been approved 
by the Executive Committee of the Washington, D. C., 
Section of the IRE and should get under way before the 
first of the year. Members in this area should contact 
Mr. Clarence Stewart of Jansky & Bailey, Inc. for fur- 
ther information. 

The Los Angeles Orange Belt sub-section has been 
initiated by A. W. Walters of Naval Ordnance Labora- 
tory at Corona, California. Members of this area should 
contact Mr. Walters for further information. 

The first regularly scheduled meeting of the Phila- 
delphia Chapter of the IRE Professional Group on An- 
tennas and Propagation was held November 17, 1954 at 
the Towne Building, University of Pennsylvania. This 
meeting was held jointly with the Philadelphia Chapter 
of the IRE Professional Group on Microwave Theory 
and Techniques. Attendance totaled twenty-two mem- 
bers and guests. 

The following officers for the Philadelphia Chapter 
were unanimously elected: Chairman: Mr. Charles Polk, 
Moore School of Electrical Engineering, University of 
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Pennsylvania; Vice-Chairman: Mr. R. A. Dibos, Philco 
Corporation; and Secretary: Mr. S. M. King, I-T-E 
Circuit Breaker Company. 

During the technical portion of the meeting two pa- 


pers were presented: 
“Shielding with Coated Glass or Thin Conductive Sur- 


faces,” by Dr. E. I. Hawthorne, University of Pennsyl- 


vania. A mathematical analysis of the transmission of 
electro-magnetic fields through thin conducting planes, 
infinite in extent, was presented. In particular, two 
sources were considered: a magnetic and an electric cur- 
rent element assumed normal to the sheet. The field 
along the axis of the source on the opposite side of the 
sheet was determined and compared with the field which 
would exist if the sheet were not present. ; 

“Theoretical Gain of Microwave Reflectors,” by D. R. 
Crosby, R.C.A. The physical behavior of reflectors from) 
the field point of view was described. Various methods of 
numerical analysis were indicated, and the results of a 
particular computation were shown. 

Both papers proved to be highly informative to the 
assembled group, and were followed by lively discus-. 
sions. After the technical papers, the chair announced | 
that the next regularly scheduled meeting would again. 
be a joint meeting of the PGAP and the PGMT&T to) 
be held at 8 P.M. on Wednesday, February 16, 1955, at) 
the University of Pennsylvania. The highlight of this: 
meeting is to be a paper by Dr. S. A. Schelkunoff en-- 
titled, “The Early History of the Theory of Dipole An-, 
tennas.” 


ANTENNAS AND PROPAGATION PEOPLE 


Mr. Kenneth A. Norton has been appointed chief of 
the Radio Propagation Engineering Division of the’ 
Boulder Laboratories of the National Bureau of Stand-. 
ards in Colorado. He was formerly chief of the Systems) 
Research Laboratory of the NBS Central Radio Propa-. 
gation Laboratory. When the CRPL moved from Wash-. 
ington, D. C. to Boulder, Mr. Norton’s laboratory be-: 
came a division. The primary responsibilities of this di-: 
vision are tropospheric propagation research and the re-. 
lated studies of frequency utilization. Mr. Norton has: 
been a member of the NBS staff since 1946. He also) 
served previously at the Bureau, in 1929 and 1930. He: 
was a consultant in the Radio Propagation Section of) 
the Office of the Chief Signal Officer from 1944 to 1946, , 
radio and tactical countermeasures analyst in the Oper-.- 
ational Research Section of the Eighth Air Force in Eng: 
land from 1943 to 1944, and Assistant Director of the: 
Operational Research Group in the Office of the Chief! 
Signal Officer from 1942 to 1943. Prior to World War II,, 
he was for a number of years on the staff of the Techni-- 
cal Information Section of the Federal Communications: 
Commission. Mr. Norton is a native of Rockwell City,, 
Iowa, and received the B.S. degree in physics from the 
University of Chicago in 1928. 

Dr. Ralph J. Slutz has been appointed Chief of the 
Radio Propagation Physics Division of the National 


1955 


Bureau of Standards Boulder Laboratories in Colorado. 
This division is a part of the Bureau’s Central Radio 
Propagation Laboratory. Dr. Slutz -will direct the Bu- 
reau’s research on the physics of radio propagation with 
particular reference to the ionosphere. Prior to joining 
the Central Radio Propagation Laboratory early in 
1954, Dr. Slutz was well known for his role in the design 
and construction of SEAC, one of the Bureau’s elec- 
tronic digital computers, and for his research on the ap- 
plication of electronic computers to the solution of 
problems in theoretical physics. From 1949 to 1953 Dr. 
Slutz was Assistant Chief of the NBS Electronic Com- 
puter Laboratory. From 1946 to 1948, he was a design 
engineer at the Institute for Advanced Study at Prince- 
ton, New Jersey, doing research on electronic computer 
design; from 1942 to 1945 he was a member of Division 
2 of the Office of Scientific Research and Development, 


investigating the effects of bomb and projectile explo-. 


sion and impact; and in 1937 and 1938 he was with Bell 
Telephone Laboratories. 

Professor Samuel Silver of the University of California 
has recently returned from a sabbatical year during 
which time he engaged in research in two European uni- 
versities. Traveling with his family on a Guggenheim 
Fellowship for Research in Diffraction, he divided his 
time between the Royal Technical Institute in Copen- 
hagen, Denmark and Cambridge University in England. 
Before returning to Berkeley he attended the eleventh 
General Assembly of URSI in the Hague, Netherlands. 
In addition to his activities in the technical programs of 
the Assembly, Professor Silver was honored by re-elec- 
tion as President of the International Commission on 
Radio Waves and Circuits (Commission VI of URSI) 
for the coming three-year period. 

Dr. E. C. Jordan, formerly professor of Electrical En- 
gineering at the University of Illinois, has been named 
Chairman of the Electrical Engineering Department of 
that University. 

Dr. V. H. Rumsey, previously head of the Antenna 
Laboratory at Ohio State University, has been ap- 
pointed professor of Electrical Engineering at the Uni- 
versity of Illinois. 

Dr. Chen To Taz has recently been appointed Associ- 
ate Professor in the Department of Electrical Engineer- 
ing at Ohio State University, where he will teach anten- 
na courses and do research in the Antenna Laboratory. 
Dr. Tai is well known for his many papers on antenna 
theory. Born in 1915 in Soochow, China, he came to the 
United States in 1943, and received his D.Sc. degree 
rom Harvard University in 1947. From 1943 to 1949 he 
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was a teaching fellow and a research associate specializ- 
ing in antennas at Harvard University. In 1949 he 
joined the Stanford Research Institute where he was 
Senior Research Physicist. Dr. Tai is a member of Sig- 
ma Xi and the American Physical Society. 

Dr A. H. Waymeck and Dr. H. G. Booker ate ona 
year’s leave of absence on a special teaching assignment 
to Cavendish Laboratory in Cambridge, England. 


CORRESPONDENCE 


In the October, 1954 issue of the PGAP Transac- 
tions, the Editor asked for opinions on several topics. 


- W. Sichak of the Federal Telecommunications Labora- 


tories has sent us his. 

“First, the Editor and the reviewers are to be com- 
mended for the decreasing time lag between receipt of 
a paper and its publication. One valid criticism of the 
PROCEEDINGS OF THE I.R.E. in previous years was that 
papers took over a year (on the average) to be published. 
With good teamwork, it should be possible to publish 
papers in less than six months, as is done by the Ameri- 
can Physical Society. 

“Papers on microwave components should not be 
published by the PGAP, but rather by the Group on 
Microwave Theory and Techniques. Papers on Radio 
Astronomy should by all means be published by PGAP, 
because the subject is closer to the interest of this Group 
than to that of any other IRE Group; the antennas and 
techniques used can be of use in other fields; and the 
data collected by radio astronomers is useful even now, 
since in the design of a modern radio relay link consid- 
eration must be given to cosmic noise sources that might 
come within the beam of the receiving antenna. 

“Good papers from our friends abroad should be en- 
couraged, since we have no monopoly on good ideas. 
Further, it might encourage more of us to read foreign 
journals and thus decrease our technical insularity. 

“The reference to special handling of student papers 
is not understood. If it means lowering the standards for 
student authors, I vote against it.” 
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Double Parabolic Cylinder Pencil-Beam Antenna’ 


R. C. SPENCER}, F. S. HOLT}, H. M. JOHANSON, anp J. SAMPSON} 


Summary—Radiation from a point source placed on the focal 
line of a parabolic cylinder is reflected in succession from this 
cylinder and from a second parabolic cylinder crossed so that its 
focal line coincides with the directrix of the first cylinder. The two 
reflections result in a parallel beam. The theory is applicable to both 
microwaves and light. The advantages of shipping the cylinders in 
the form of flat sheets and the possibilities of independent control of 
horizontal and vertical beamwidths and shapes are pointed out. Ex- 
perimental models have been built and tested. 


INTRODUCTION 


N CONTRAST to the usual microwave pencil- 

beam antenna formed by allowing energy from a 

point source to reflect once from a paraboloid of 
revolution, the device to be discussed makes use of suc- 
cessive reflections from two parabolic cylinders with ele- 
ments at right angles, each cylinder collimating the beam 
in one plane. Functionally, the point source and first 
parabolic cylinder could be replaced by a parabolic 
pillbox! or any other equiphase line source. 

The idea of using crossed cylindrical lenses has long 
been known in optics’, but the application to pencil 
beams appears to be new. More recently, the use of 
crossed cylindrical reflectors has been introduced in 
connection with X-ray microscopes**, where the phenom- 
enon of total reflection of X-rays at grazing incidence 
lends itself to the use of curved metal reflectors. 


* Original manuscript received by the PGAP, June 21, 1954; 
revised manuscript received, October 11, 1954. The geometric design 
was discussed under the same title, by RG Spencer, in Air Force 
Cambridge Research Center Rep. No. E5084; April, 1952. 

+ Air Force Cambridge Research Center, Cambridge, Mass. 

t Now in the U. S. Navy. 

1S. Silver, “Microwave antenna theory and design,” M.I.T. 
Rad. Lab. Ser., vol. 12, pp. 151, 494; 1949. 

2 J. P. C. Southall, “Mirrors, Prisms and Lenses,” The Macmillan 
Company, New York, N. Y., ch. 1X; 1933. 

EM ee Kirkpatrick and A. V, Baez, “Formation of optical images by 
X-rays,” Jour. Opt. Soc. Amer., vol. 38, pp. 766-774; 1948. 

aie Ehrenberg, “The Production of converging beams by total 
reflection,” Jour. Opt. Soc. Amer., vol. 39, pp. 741-746; 1949, 


Kirkpatrick and Baez? discuss the use of two crossed! 
elliptic cylinder mirrors to form a real image of a point. 
and the use of two parallel parabolic reflectors to form a} 
real line image of a point. In references 3 and 4, the re-- 
flectors are approximated by sections of circular cylind-. 
ers. In 4, the elliptic cylinder is approximated by a bent: 
optical flat. 


GEOMETRY 


Fig. 1 (opposite) shows a perspective drawing of the: 
two parabolic cylinders. The first cylindrical surface Su 
is given by the equation 


2 = Afx = 2x, (1)) 


where f is the focal length and /=2f is the semilatus: 
rectum. The point F with co-ordinates (f, 0, 0) lies on) 
the focal line of S;. The line DD’ is the image of point F' 
in S$); 1.e., any ray emitted from F and incident on S; at: 
P, is reflected as though it originated on DD’. | 

The formula for the second cylindrical surface S: is: 
given by the equation 


(a +f)? = 4f@ +f). (2) 


This surface has been so oriented that it meets the fol-. 
lowing conditions: (a) its focal line coincides with DD’, 
(b) it contains the point F, and (c) any ray emitted 
radially from DD’ is reflected parallel to the positive g 
axis. Under these conditions, S,; has the same focal! 
seh f as S;, and the point V with co-ordinates (—f, 

0, —f) lies on the vertex line of Sy. 

It is evident, therefore, that any ray FPP, emitted 
from a point source located at F will, after successive 
reflections from S, and 53, be directed parallel to the 
positive z axis. If either (b) or (c) on S_ were relaxed, 
then other combinations of relative focal lengths and 
orientations would be possible. 
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_ Fig. 1—Orthogonal curvilinear co-ordinate systems on surfaces 
S; and So. 


CALCULATION OF THE CO-ORDINATES OF 
COMMON INTERSECTION 


The x, y, and zg co-ordinates of the intersection of the 
‘wo parabolic cylinders with 2f=/=1 were calculated 
ising (1) and (2). The curve of intersection may be 
vritten in the following parametric form: 


- Spaces \ 
5 ae 


es | e 


8 2 


a 
| 


Table I contains numerical values of x, y, and z as 
omputed from (3). The range on ¢ (that is, y) is from 
) to 2 in increments of 0.1. 


METHOD OF DEVELOPING CYLINDERS 
ONTO FLAT SHEETS 


The curve of intersection of the two surfaces S; and 
$2 is a three-dimensional space curve that we shall de- 
iote as C (see Fig. 1). If the surfaces S; and S» are sepa- 
ately developed into planes, the intersection curve in 
ach case becomes a plane curve; however, these two 
lane curves are not identical. In the construction of 
he double cylindrical reflector, it is desirable to know 
he co-ordinates of these plane curves in order that they 
nay be laid out on flat sheets before the sheets are bent 
nto parabolic cylinder shape. 
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Curves of z=constant in the surface S are all iden- 
tical parabolas and S; is defined to be the are length 
along any of these parabolas as measured from the 
2 axis (see Fig. 1). Similarly, the curves of y=constant 
in the surface S) are all identical parabolas and sz is de- 
fined to be the arc length along any of these parabolas 
as measured from the vertex line HH’. In Table I we al- 
ready have the x, y, and z co-ordinates of the intersec- 
tion curve C. In order to obtain the s,, z co-ordinates of 
C as measured in the surface |S}, it is necessary to deter- 
mine s; as a function of either z or y; that is, the arc 
length of a parabola as a function of its co-ordinates. 
Similarly, in order to obtain the 5», y co-ordinates of C 
as measured in the surface Ss, it is necessary to deter- 
mine s2 as a function of either x or z. Tables of parabolic 
arc length vs co-ordinates were computed by the Center 
of Analysis, Massachusetts Institute of Technology.® 
The last two columns of Table I were obtained from 
these parabolic arc length tables. 


TABLE I 


COORDINATES OF THE CURVE OF INTERSECTION OF 
THE TWO PARABOLIC CYLINDERS 


Point t 
NG: x y Zz S1 52 
0 0.0 0.000 0.0 |—0.3750 0.0000 0.5201 
1 ail .005 eS 3 725 .1002 SAS 
2 oP .020 2 |— .3648 2013 5426 
3 ae) .045 SOP 10515 3044 .5709 
4 A .080 4 |— .3318 .4104 .6110 
5 0.5 ORS 0.5 |—0.3047 0.5201 0.6636 
6 6 .180 .6 |— .2688 6343 .7293 
7 si 245 = SYS) .7536 8091 
8 8 .320 8 |— .1638 .8786 9043 
9 9 405 9 |— ,0905 1.0098 1.0166 
10 1.0 0.500 1.0 0.0000 1.1478 1.1478 
11 (hoa .605 1.41 .1105 1.2928 1.3003 
12 hee .720 2 2442 1.4452 1.4766 
13 123 .845 dnd 4045 1.6053 1.6799 
14 1.4 .980 1.4 .5952 1.7733 1.9136 
15 15 i SWS) 125 0.8203 1.9495 2.1814 
16 1G 1.280 1.6 1.0842 2.1339 2.4874 
tk se 1.445 il oa 1.3915 2.3269 2.8363 
18 thet} 1.620 1.8 1.7472 2.5284 Se 2027 
19 1.9 1.805 1.9 2.1565 2.7387 3.6819 
20 270 2.000 2.0 2.6250 2.9579 4.1893 


The straight-line elements of the cylinder S; along 
which the g co-ordinate is measured, together with the 
parabolas along which s; is measured, form an orthogo- 
nal curvilinear co-ordinate system on the surface Sj. 
When the cylinder is developed into a plane, this co- 
ordinate system becomes a rectangular system in the 
plane. The si, z co-ordinates of the intersection curve C 
are given in Table I. Fig. 2(a), next page, shows a plot 
of this curve on the developed surface 5S}. Similarly, the 


5 R. C. Spencer and G. E. Reynolds, “A Table of Normalized 
Parabolic Coordinates and Arc Lengths,” Air Force Cambridge Re- 
search Center Report No. E4083; July, 1951. 


Fig. 2—(a) Curve of intersection developed on surface .S;. 
(b) Curve of intersection developed on surface S». 


cylinder Sz contains an 52, y orthogonal curvilinear co- 
ordinate system that becomes a rectangular system 
when 2 is developed into a plane. The s2, y co-ordinates 
of the curve of intersection are given in Table I, and a 
plot of this curve on the developed surface S_ is shown 
in Fig. 2 (b). Note that certain points along the inter- 
section curves, as shown in Figs. 2 (a) and (b), are num- 
bered according to the first column of Table I. When 
the developed surfaces S, and S; are returned to their 
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parabolic cylindrical form and fitted together, identi- | 
cally numbered points will coincide. 


INVESTIGATION OF FEED ORIENTATION FOR 
ZERO CROSS-POLARIZATION 


Assume an electric dipole located at the focal point F 
(see Fig. 3) and polarized in the direction of the unit 
vector 


(4) | 


where i, ji; and & are unit vectors in the x, y, and g di- | 
rections respectively. The directions of the rays leaving 
the source are determined by the angles ¢ and 6, meas- 
ured at the point F. The angle ¢ is measured between 
the xy plane and the ray, and the angle @ is measured — 
between the negative x axis and the projection of the | 
ray on the xy plane. The electric field vector along any | 
ray is initially assumed to be directed normal to the ray | 
and to lie in the plane of the @ vector and the ray. Satis- | 
faction of the boundary condition that the tangential _ 
component of the electric field vanish on each reflecting | 
surface leads to the following expression for the electric | 
field on the aperture plane of the double parabolic — 
cylinder antenna: 


1 


A =ia+ jb ke, 


ua 


: {i[(a cos 6 — b sin 8) sin@ +c cos ¢| 
sin 


+ j[a sin @+ b cos @|}. 


(S) 


Fig. 3—Feed orientation and polarization. - | 


If there is to be no cross-polarization, that is, if E is to | 
have a constant direction, the components of E must be. 
proportional, as follows: | 


K [(a cos 6 — 6 sin @) sin ¢ + ¢ cos $| 


=asin@+ 6 cos @. | 


The independence of the variables 6 and ¢ and the linear 
independence of the functions sin @ and cos @ lead to the 
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following result: 
a=b=K=0. 


Since es =1, we have c=1; and (4) shows that this re- 
quires the source to be polarized parallel to the z-axis. 


ILLUMINATION OVER THE APERTURE 


An exact ray-tracing analysis of the antenna leads to 
a geometric optics approximation of the aperture illum- 
ination. Assuming an isotropic source at the feed posi- 
tion, the illumination over the aperture is found to be 
proportional to a function separable in the normalized 
co-ordinates x/f and y/f as follows: 


P= kF(x/1)G(y/f), (6) 
1 x / f 
F(«/f) = 24.63 ra i (7) 
4 
G ee eee 8 
CU Ea eaeraare (8) 


Fig. 4 shows plots of these normalized functions. 


Sele 


1.0 


yt 
(b) 
Fig. 4—Aperture illumination. (a) « direction. (b) y direction. 


EXPERIMENTAL MODEL DESIGN 


An experimental model of the double parabolic 
cylinder antenna with aperture 1 foot by 1 foot, f=6 
peches was designed and built. Full use was made of the 


Spencer, Holt, Johanson, Sampson: Double Parabolic Cylinder Pencil-Beam Antenna i 


table computed by Spencer and Reynolds? in calculat- 
ing the shapes of the flat sheets that were later bent to 
form the parabolic cylinder surfaces. In Fig. 5 (a) the 
reflector surfaces and the component parts of the sup- 
porting framework are shown laid out prior to assembly. 
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Fig. 5—Experimental model. (a) Disassembled. (b) Assembled. 


The assembled antenna is shown in Fig. 5 (b). The me- 
chanical design in this particular model features com- 
pact packaging in the disassembled state as well as 
simple rapid assembly and take-down procedures 
through use of dowels and slide snap fasteners. 

The amplitude-tapering effects of the system as shown 
in Figs. 4 (a) and (b) were taken into account in the de- 
sign of a K-band (A=1.25 cm) horn feed that theoreti- 
cally produced at least a 10 db illumination taper over 
the final aperture in both the x and y directions. It is 
evident from Fig. 4 (a) that if the primary pattern of 
the feed were symmetric in ¢ (see Fig. 3), then the final 
aperture illumination would be asymmetrical in the x 
direction. This asymmetry in illumination introduced 
purely by the geometry of the system was counteracted 
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(b) 
Fig. 6—Experimental patterns taken with the K-band model. 


to some extent by orienting the feed horn so that its 
primary pattern was asymmetrical in @. 

Experimental patterns taken with the K-band model 
are shown in Fig. 6 (a) and (b). It is felt that the asym- 
metry of the side lobes in the E-plane pattern is the re- 
sult of interference by the reflector S; with the final col- 
umnated radiation from reflector S: (see Fig. 5 (b)). 


SHAPED BEAM DESIGN (COSECANT SQUARED) 


After reflection from the first parabolic cylinder Si, 
all radiation striking the second parabolic cylinder S_ ap- 
pears to originate from a line source located along the 
line DD’ (see Fig. 1, page 5). Hence it is necessary to re- 
design only the cylinder S: in order to produce a shaped 
fan beam. 

Using the method suggested by L. J. Chu,® the design 
of the second cylindrical surface S; of the experimental 

6 R. C. Spencer, “Synthesis of Microwave Diffraction Patterns 


with Application to csc? @ Patterns,” Mass. Inst. Tech. Rad. Lab. 
Report No. 272(54-24), pp. 5-7; June 23, 1943. 


Fig. 7—Experimental elevation and azimuth patterns. 


model was modified to produce a csc? 6 fan beam over 
the range 2 to 20 degrees. Experimental elevation and 
azimuth patterns are shown in Fig. 7 (a) and 7 (b). 


CONCLUSION 


The double parabolic cylinder antenna has several ad- 
vantages over the conventional paraboloid. Each of its 
surfaces, being cylindrical, can be stamped from flat 
sheets. The system itself can be designed so that it can 
be easily dismantled and stored or transported in a 
compact flat package. There is an absence of cross- 
polarization. Since the first reflector effectively focuses 
in one plane and the second in the orthogonal plane, a 
degree of independent control can be obtained over the 
beam shapes in both planes by modifying the appropri- 
ate reflecting surfaces. In particular, the second surface 
S. can be modified to obtain csc? 6 patterns. 
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An Atmospherics Analyzer* 
PERRIN F. SMITHt 


Summary—The enormously powerful radiation from lightning 
discharges has been used as a tool for the study of the ionosphere at 
very low frequencies. Previous workers have calculated the height 
of the ionosphere and the distance to the radiation source by appli- 
cation of the ray theory of multiple ionospheric reflections. Time 
measurements for use in the height and distance formulas were ob- 
tained by direct measurement of recorded oscillograms. Time reso- 
lution was limited and the calculations were extremely tedious. 

This paper describes a complete system which records the re- 
ceived atmospheric wave form and simultaneously derives a plot of 
the time intervals between peaks of the wave form vs time. The time- 
interval plot greatly facilitates data reduction by means of a simple 
curve-fitting technique. Accuracy of measurement of height and dis- 
tance is limited chiefly by the validity of the multiple reflection the- 
ory and the assumption of a simple layer structure in the ionosphere. 
The time differences are measured more accurately than is ordinar- 
ily possible with direct measurement on the waveform. The system 
automatically derives frequency-time plots for whistling atmospher- 
ics and relates the whistler to preceding impulsive atmospherics. 
Consumption of photographic film is low. The system can be used to 
visually monitor tape-recorded atmospherics. 


INTRODUCTION 


HE ENORMOUSLY powerful radiation from 

lightning discharges is potentially a powerful tool 

for the study of the ionosphere at very low fre- 
quencies. Analysis of certain types of recorded wave- 
forms will indicate the average height of the ionosphere 
over the transmission path and the distance to the dis- 
charge. 

This paper describes a recording system for atmos- 
pherics which presents data from the received wave- 
forms in a partially reduced form, greatly facilitating 
the determination of the height of the ionosphere and 
the distance to the discharge. The time resolution of the 
system is increased several times over the resolution 
obtainable with a single-sweep display. The ease of data 
reduction with the analyzer described should make prac- 
tical a statistical study of a large number of waveforms. 


IONOSPHERIC PROPAGATION OF ATMOSPHERICS 


For the commonly occurring case of fairly small losses 
in the ionosphere, simple ray theory is sufficient to ex- 
plain the important features of the observed wave- 
forms.! The single pulse of radiation from the lightning 
discharge is received at a distant point a number of 
times, each time corresponding to a particular number 
of hops between earth and ionosphere. The ionosphere 
is assumed to be a smooth perfectly reflecting layer. 


* Original manuscript received by the PGAP, February 19, 1954; 
revised manuscript received July 20, 1954. t 

This paper was presented at the joint meeting of the Inter- 
national Scientific Radio Union, USA Natl. Com. and Can. Natl. 
Com., and the IRE PGAP in Ottawa, Can.; October 6, 1953. 

+ Associate Engineer, International Business Machines Corp., 
Engineering Laboratory, San Jose, Calif. ! af 

i K, G. Budden, “The propagation of a radio-atmospheric, The 
Philosophical Magazine, vol. 42, pp. 1-19; 1951. 


Budden!” has proposed a general theory to explain the 
propagation of atmospherics initiated by lightning dis- 
charges. 

The transmission time of the mth-order pulse may be 
shown from the geometry to be 


2n D 
ia i y/ 2R(R + h) (1 — cos ) + h*,. ) 
C 2Rn 


where »=number of reflections at the ionosphere, 
t,=time of arrival of the mth impulse, c=velocity of 
light, h=virtual height of reflection at the ionosphere, 
R=radius of earth, and D=great-circle distance to 
radiation source. 

If cos D/2Rn is approximated by the first two terms 
of the series, we have very closely 


1 h 
ie =4/o(1 + =) + 4n7h?. (2) 
C R 


When the order of reflection is known, the following 
expressions for the height of the ionosphere and the dis- 
tance to the source can be derived from (2): 


T,T2(T TE 
ies c / 1 2( ene 2) (3) 
ON AG Cte) tO mei heh AE 
C nee p?) an ae sp 
2cT, 
D? = (4) 
ieee 
R 


where , g, and r are pulse orders and 
10; = be = bp 
dis, = ie = tae 


The first two or three pulses arriving at a receiver from 
a long distance tend to merge, so that the first clearly 
defined pulses received are usually of the second or third 
order.’ 

Others*® have made use of (3) and (4) in determining 
h and D by directly measuring the recorded waveform. 
Accuracy of measurement is very important, as is evi- 
dent from (3), because the height measurement depends 
on a difference of differences. Order must be determined 


2K. G. Budden, “The propagation of a radio-atmospheric II,” 
The Philosophical Magazine, vol. 43, pp. 1179-1200; 1952. 

3 For the adverse but typical case of D=2000 km and h=90 km, 
the exact result for 4, differs from the approximation found from (2) 
by about 1.5 per cent. The approximation improves rapidly with 
increasing ”. The chief limitations on accuracy are the validity of 
the simple ray theory and experimental error. 

4P. G. F. Caton and E. T. Pierce, “The waveforms of atmos- 
pherics,” The Philosophical Magazine, vol. 43, pp. 393-409; 1952. 

5 W. J. Kessler and S. E. Smith, “An atmospheric waveform 
receiver,” Proc. I.R.E., vol. 39, p. 676; June, 1951. 
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by a trial-and-error procedure and h and D should be 
found for several groups of pulses to obtain average 
results. If many waveforms are to be analyzed, the data 
reduction problem becomes difficult. 

Curves of the time separation between pulses of the 
received waveform vs time may be prepared from (2) 
for a range of values of # and D. Each curve will have a 
characteristic shape. The atmospherics analyzer gen- 
erates a dot pattern on an oscilloscope from each re- 
ceived wave form which is essentially the At vs ¢ plot 
from (2). By curve-fitting on a photographic record of 
the generated dot pattern it is possible to select very 
rapidly the “best” height and distance curve from a 
family based on (2). 


TEKTRONIX 


511 


SWEEP 
OUTPUT 


OSCILLOSCOPE 


INTER VALOMETER 


Fig. 1—Block diagram of recording system. 


A block diagram of the recording system is shown in 
Fig. 1. The loop antenna is mounted in a manner to 
permit rotation to null out power line hum or interfering 
stations. The wide-band amplifier is of conventional 
design and has a pass band from approximately 1 kc to 
20 kc. The limited pass band was necessitated by inter- 
ference from nearby low-frequency transmitters. The 
resolution of the system was not significantly increased 
by including frequency components above 20 kc; in fact, 
an upper frequency limit of 10 ke was found adequate. 

The auxiliary trigger device is provided to assure that 
the sweep will be initiated at the beginning of a received 
waveform regardless of the sign of the first half-wave. 

The sweep output of the synchroscope is direct- 
coupled to the y-axis input of the direct-coupled oscillo- 
scope (a DuMont 304H). The two oscilloscopes then 
sweep together. 

The intervalometer is essentially a linear-sweep gen- 
erator which re-initiates its sweep upon the occurrence 
of each input pulse. The time separation between pulses 
is thus measured by the excursion of the corresponding 
sweep. The reference point for time measurement on 
each pulse is taken as the positive peak of the pulse, 
identified as the point of zero first derivative. 

The operation of the intervalometer can be under- 
stood by reference to Fig. 2, which shows the separate 
waveform operations in detail. 
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The input pulse is amplified to a convenient level, 
then base-clipped at a positive voltage sufficiently high 
to eliminate low amplitude components of the pulse 
and interference. The positive portion of the clipped 
waveform is differentiated, peak-clipped, and differ- 
entiated again to obtain a narrow pulse at the point of 
zero derivative of the input pulse. This derived pulse 
triggers a delay multivibrator, which triggers a second 
multivibrator which delivers a large positive pulse to a 
sweep circuit. The sweep circuit is simply a power triode, 
biased to cutoff, which has a condenser between plate 
and ground. An input pulse to the sweep tube causes it 
to conduct heavily and discharge the sweep condenser 
to zero, assisted by stray inductance in the circuit. 
The circuit must operate on irregularly occurring in- 
puts, hence must be capable of discharging the sweep 
condenser and its output-coupling condenser from the 
power supply voltage to zero in a time less than the 
minimum expected interval (about 30 wsec) between 
input pulses. 

The positive pulse from the delay multivibrator is 
used as a z-axis intensifying signal to the direct-coupled 
oscilloscope. A short bar of light then appears on the 
oscilloscope screen at the time of the peak of each input 
pulse. At the end of the short intensification period, the 
sweep circuit discharges to zero. 
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Fig. 2—Waveform operations of intervalometer. 


The presentations of both oscilloscopes are recorded 
ona single 35-mm film strip by means of the camera 
and mirror system shown in Fig. 1. The shutterless 
camera is driven at approximately 2 inch per second. 
Fig. 3 (opposite) illustrates the nature of the presenta- 
tion. The received waveforms are recorded in the lower 
half of the picture, time running downward the output 
of the intervalometer, in the upper half of the picture, 
increasing time difference running from the center up- 
ward. The long-time sweep runs right to left and is 
approximately 9 msec in duration. Film motion is op- 
posite to the sweep motion. Note that the origin of the 
atmospheric waveform and the origin of the locus of 


points derived from the waveform are displaced from 


each other on the record. A rear-surfaced mirror was 
used in making these pictures, hence the faint “ghost” 
in the waveform record which must be disregarded. 
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Fig. 3—Typical good atmospheric waveform and derived locus. 


In studying atmospherics, it has been found con- 
venient and economical to record the atmospherics on 
magnetic tape and to operate the system with the tape 
record. Visual monitoring of the intervalometer pres- 
entation is used to decide what portions of the tape 
record are of sufficient interest to photograph. 


SAMPLE RECORDS AND DISCUSSION 


Fig. 3 is a good example of a locus resulting from an 
atmospheric waveform which is well explained by the 
simple theory of (2). One curve of time separation be- 
tween pulses vs total time may be fitted to the locus of 
points from a family of curves previously computed 
from (2). There will be one best fit for a given height 
and distance.® The locus and waveform on the right side 
of Fig. 4 illustrates the limitations of the assumption of 
a simple layer structure in the ionosphere. 

The atmospheric waveform indicated on the right 
side of Fig. 4 is of the type classified by Caton and 
Pierce* as “quasi-sinusoidal.” Note how a locus is de- 
fined in the presence of many waveforms occurring at 
nearly the same time. Note also that the locus in this 
case departs from the smooth shape predicted by (1). 
The reason for the slightly erratic pattern is apparent 
from the skewed shape of the pulses in the wave form 
itself. The act of fitting a curve to the locus accomplishes 
essentially the same result as finding h and D from (3) 
and (4) for several intervals measured on a waveform 
and averaging the results. The erratic points at the end 
of the intervalometer locus correspond to the distorted 
pulses at the end of the waveform. The best fit to this 
rather poor locus is h=90 km and D=1,500 km. 


6 Note that # and D are not found by direct application of (3) 
and (4) but by identifying a particular change of time interval with 
time. Constant errors in measuring the time intervals will shift the 
locus slightly, producing errors in curve matching roughly propor- 
tional to the shift. Application of (3) and (4) involves the difference 
of differences between intervals, and small errors in reading individual 
intervals can produce large errors in h and D. 
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Fig. 4—Atmospheric waveforms followed by whistler. 


Curve-fitting by eye has the advantage of much 
greater speed and also permits one to reject points on a 
locus which apparently do not belong to the locus or to 
give them a weighted value. About seven points are 
necessary to identify a locus with a curve. The accuracy 
is chiefly limited by the quality of the locus obtained. 


Fig. 5—Multiple atmospherics. 


Good loci should give height to +3 km and distance to 
+50-200 km, depending upon the distance. Poor loci 
should give two or three times this error. The group of 
multiple atmospherics on the left side of Fig. 5 is prob- 
ably of the same origin as in Fig. 4, observed at the 
same time. It is difficult to fit curves accurately to loci 
from rapid multiples at low film speeds. 

The atmospheric discussed above was followed ap- 
proximately one second later by a large whistler’ (see 
left side of Fig. 5). This whistler was of sufficient ampli- 


7A whistler is a relatively rare atmospheric which is usually 
heard as a musical tone descending from approximately 5 ke to 
about 1.5 ke in a period between 0.25 second and 1.5 seconds. The 
musical tone often has a swishing sound. 
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tude to trigger the synchroscope almost continually. 
Operation of the intervalometer during the dead time 
between synchroscope sweeps caused the bright bands 
to be recorded as in Fig. 4 with the time base provided 
by film motion. Note that a second weaker band is re- 
corded which diminishes in frequency more rapidly than 
the brighter band. This effect has been observed by 
Potter,’ but not while simultaneously recording the 
atmospheric preceding the whistler. The author has ob- 
served the tendency to form discrete bands which 
change frequency at different rates in most of the larger 
whistlers studied. The second band appears to diminish 
in frequency roughly twice as fast as the main bright 
band. Many whistlers appear to have maximum ampli- 
tude at the asymptotic limiting frequency of the large 
atmospheric that often precedes the whistler. 

An “integrated locus” characteristic of an electrical 
storm giving rise to atmospherics may be obtained by 
making exposures of several seconds without film mo- 
tion. With good conditions, the superposition of the loci 
from many individual atmospherics will form a single 
clearly defined locus which is characteristic of the layer 
height and the distance to the storm. 


CONCLUSION 


The recording system described should be a useful 
tool for the study of the ionosphere at very low fre- 


8 R. K. Potter, “Analysis of audio frequency atmospherics,” 
Proc. I.R.E., vol. 39, pp. 1067-1069; September, 1951. 


quencies by the use of atmospherics. With the addition 
of a direction-finding function, sources of lightning 
activity might be located and tracked. The ease of ob- 
taining height and distance data makes practical the 
statistical study of large numbers of waveforms. Great 
economy is obtained by first recording the atmospherics 
on magnetic tape. Following that the interesting parts of 
the record are visually selected for transcription to film. 


Two improvements of the recording system for atmos-_ 
pherics described in this paper are suggested. The first ) 
is to produce points on the locus from both positive and_ 


negative half-cycles of the atmospheric waveform. The 
second improvement is to make the base-clip level for a 
given pulsedepend on theamplitude ofthe preceding pulse. 
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A Single-Control Tuning Circuit for 
Electrically Small Antennas” 
R. E. WEBSTER} 


Summary—lIt is often desirable to tune electrically small anten- 
nas at their feed terminals. An ‘‘L” matching circuit for this purpose 
is considered. It is shown that one fixed element and one variable 
tuning element are adequate for certain small antennas. The manner 
in which losses affect the circuit effiiciency and the single-control 
tuning feature is discussed. Application of the technique to a tunable 
cavity antenna is noted. 


INTRODUCTION 
Peiienae sonca small antennas may be con- 


sidered as nonresonant high-Q two-terminal net- 
works. When such antennas are resonated with 
an ideal tuning element, the circuit half-power percent- 
age bandwidth is given by 
Ee Original manuscript received by the PGAP, June 21, 1954. 
This article is based on work done under a contract between the 
Ohio State Univ. Res. Foundation and the Signal Corps Eng. Labs., 
Ft. Monmouth, N. J. 


rit Research Associate, Antenna Lab., Ohio State Univ., Columbus, 
io. 


1 Re Gs. : 
HPBW = —= = 
Q |X| |B| 
where the antenna impedance is 
1 
Za = hee are ee ar >; ; 
Ye Got7Ba 


It is well appreciated that variable tuning is necessary | 
if an impedance match is desired over a frequency 


bandwidth larger than a few per cent. A minimum of 


two tuning elements is required to match Z, to a fixed | 
real impedance, and generally both elements must be | 
variable. The following comments concern a simple | 


“L” matching circuit which, when used with impedances 


ime’ Montgomery, R. H. Dicke, and E. N. Purcell, “Principles 
of Microwave Circuits,” McGraw-Hill Book Co., Inc., New York, 
N. Y., Sec. 4-13; 1948. 


| 
} 
} 
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ybtained with small antennas, can consist of one fixed 
ind one variable tuning element. The required antenna 
s one for which Ra= Kif? or Ga = Kof? (the K’s are arbi- 
rary constants and f denotes frequency). This fre- 


juency dependence is realized in certain small dipoles, 
oops, and slots. 


DEFINITIONS 


The electrically small antennas considered are (1) 
lipoles less than 0.20 wavelength long, (2) loops less 
han 0.10 wavelength in diameter, and (3) slots less 
han 0.40 wavelength long. The frequency dependence 
yf Ra or Ga mentioned above is satisfied approximately 
or these dipoles and loops, and for some slot configura- 
ions of these dimensions.?3 

The term “percentage bandwidth” will be used 
hroughout to indicate 


ic oe oa 
Q |x| |B 


when the antenna is resonated, the tuning elements be- 
ng fixed. This is the fixed tuned bandwidth of the an- 
enna and is not to be confused with the “tunable band.” 

The term “tunable band” is defined as the frequency 
‘ange over which the antenna reflection coefficient is re- 
luced to a specified value by proper adjustment of one 
uning element. 


The “tuning efficiency” 


is defined as follows: 


yercent power delivered to antenna 
tuning efficiency — 


: ee 100 
power input to tuning circuit 


The tuning efficiency is maximum when the tuning cir- 
suit losses are minimized. Since the losses in lumped 
sircuit elements are proportional to the energy storage, 
he most efficient circuit will have the least stored 
‘nergy. Thus a highly capacitive antenna will be tuned 
nost efficiently by the minimum inductive reactance 
equired for resonance. The circuits discussed here will 
ichieve this maximum tuning efficiency when they are 
-omposed of lumped elements. 


SINGLE-CONTROL TUNING CIRCUITS 
Basic Circuits 


The approach will be to combine the Af? frequency 
lependence of R, or G, with the reactance of one vari- 
ible tuning element in such a manner that the second 
equired tuning reactance will have the form X2=Kof 
x X2=—K,/f (a lumped inductance or capacitance, 
espectively). Thus the second tuning element may re- 
nain fixed. 


2H. A. Wheeler, “Fundamental limitations of small antennas,” 
Roc. I.R.E., vol. 35, pp. 1479-1484; December, 1947. 

3M. H. Cohen, T. H. Crowley, and C. A. Levis, “The aperture 
dmittance of a rectangular waveguide radiating into half space,” 
roject Rep. 339-22, Antenna Lab., The Ohio State Univ. Res. Found., 
repared under Contract W33(038)ac21114 with Air Res. and Dev. 
‘ommand, Wright Air Dev. Center, Wright-Patterson Air Force 
sase, Ohio; November 14, 1951. 
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The circuit of Fig. 1 is an equivalent circuit proposed 
for the Tunable Cavity Antenna developed at this lab- 
oratory and described at the 1953 I.R.E. National Con- 
ference on Airborne Electronics. Constructional and 


Fig. 1 


operational details of this antenna are given in a tech- 
nical report.> The assumptions for the circuit are that 
Ga = Kof? and Ga|Ba|. The expression for Zin is 


1 
Zin = 1X =F 
COUN Teac Hoses 
Ga 
Ga aie (By ae Bale 
(By + Ba) | 
+7] Xe -— : 1 
ifs Gi? + (Bi + Ba)? Y 
Let 
Ga 
= Rew (2) 
GED 2B.) 4 ; 
and 
B 18, 
g=e en) = (0), so that Z;, 1s a constant 
Gt ae ea) F 
real impedance. (3) 
Then, from (2), 
Ga 
(By + Ba) = aE x al Ga (4) 


Substituting (4) in (3), 


XQe sls 


4/ sR (5) 
Ge ah 


Eq. (3) is accomplished by properly adjusting B,, and 
(5) results. With the additional condition Gin =1/Rin> 
G., (5) becomes 


Ris. Ine 1 
LOST 4/ = + Sao errr) (6) 
Ga Kof? 2nfC 


4D. R. McCoy, “Small tunable cavity antenna for nominally cir- 
cular polarization,” Airborne Electronics; 1953. Abstract of papers de- 
livered at the 1953 I.R.E. Nat. Conf. on Airborne Electronics, 
Dayton, Ohio; May 11-13, 1953. ‘ 

5 R. W. Masters and D. R. McCoy, “A tunable cavity antenna,” 
Tech. Rep. 486-13 (Confidential) Antenna Lab., The Ohio State 
Univ. Res. Found., prepared under Contract AF 18(600)-85, Air 
Res. and Dev. Command, Wright Air Dev. Center, Wright-Patterson 
Air Force Base, Ohio; November 5, 1953. 
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where 
Ko 


en "a Re 


Thus Xe, when taken to be negative, can be realized as 
a fixed capacitor. Then, in (3), (Bi+B.) must also be 
negative. If Rin>1 (a typical value may be 50Q), and 
since | Ba| >G,, it is seen in (2) that B, must be of op- 
posite sign from B, and of somewhat different magni- 
tude. Thus, if B, is negative, B, must be positive and 
slightly smaller than B,. The result is a single-control 
tuning circuit for use with small loops and slots subject 
to the following conditions: 


(a) Ga = Kof? 

(b) GuacGn 

(c) (Bi + Ba) = — V/G.Gmn — G,? 
| eens 

i) dest ty Vela 


B, isa variable capacitor, X2 is a fixed capacitor. 

The same argument can be applied to the dual of this 
circuit (see Fig. 2). In this case the conditions for a 
constant real input admittance (Yjn=Gin) are: 


(a) Re = IGP 

(b) Ree Re 

(c) (Xi +X.) = — VRiRin — R? 
cae 

@) rte / he 


The antenna impedance is assumed to be capacitive 
with | Ro| VX) The resulting circuit is a fixed in- 
ductor Byand a variable inductor X suitable for tuning 
short dipoles. 
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first the tunable band. If loss factors Re and Gi are 
assigned to the elements of Fig. 1, modified (5) is 


Rin a R ) 
tes \ eR 

(Ga AE G) 
If G\<«G, and R:&Rin, X2 will retain the form of a fixed 
capacitor with the desired //f variation. As frequency is 
decreased, however, G.—0 as f?, and may actually be- 


Ka) a 
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(7) | 


| 


come smaller than G,, thus producing a frequency de- | 
pendence of X>» different from //f. A similar difficulty | 


could be produced by a large Ro. A fixed capacitor may 
still be used for X2, although the tunable band will be 
somewhat reduced. Thus the low-frequency end of the 
tunable band can be limited by circuit losses while the 


high-frequency limit is specified by the conditions (b). | 


If the losses are appreciable over the entire theoretical 


tunable range, of course, the actual tunable band be-— 
comes very small. These effects apply equally to the | 


dual circuit of Fig. 2, where modified (5) becomes 


(Gin re G2) 


1B = ee = (G5, = Ge) 
: (RR \ Z 


(mhos) 


REQUIRED SUSCEPTANCE 


LENGTH OF STUB (WAVELENGTHS) 


Fie. 3 


Measurements made in the 250- to 600-mc range with 


(8) 


coaxial-line tuning elements connected to a short verti- | 
cal stub protruding from a horizontal ground screen 
demonstrated these effects. The circuit of Fig. 2 was | 


Lossless tuning elements have been assumed. On this 
basis the tunable band is limited only by the percentage 
bandwidth desired at low frequencies and by the degen- 
eration of conditions (a) and (b) at high frequencies 
(when the antennas become electrically large). In prac- 
tice, however, circuit losses will considerably affect the 
tunable band as well as the tuning efficiency of the sys- 
tem. 


Effect of Losses 


Circuit losses should reduce the tunable band and the 
tuning efficiency, especially since the antennas of inter- 
est present very small radiation resistances. Consider 


simulated. Elements X; and Bz were adjusted at each 


frequency to obtain a match to a 50-ohm slotted line. | 


The tuned condition was determined by observing the | 


standing wave on the slotted line. A directional coupler — 


between the signal source and slotted line facilitated 
the tuning adjustments. Impedances of the tuning ele- 
ments were then computed from measured positions of 
their shorting plungers. Fig. 3 shows a plot of the meas- 


ured By as a function of the electrical length of the stub. 


R, was the major loss resistance. The curve is of the 


form predicted by (8) when R,; and R, dominate at the. 


low and high frequencies, respectively. The ideal varia-_ 


tion of Bz, assuming lossless elements, is also shown in 
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ig. 3; it corresponds in frequency dependence to a 
xed lumped inductor. The measured By will approach 
he ideal variation as the circuit losses approach zero. 
n the presence of losses, however, it is apparent that a 
xed inductance can be used to simulate the measured 
3o,0ver a reduced tunable band. The resulting reflection 
oefficient, | p| , ata particular frequency is 


AB | 
| V(2 Vo)? + (ABz)? 


there AB, is the difference between measured By [or Be 
s calculated by (8)], and B, of the fixed inductor. The 
orresponding expression for the circuit of Fig. 1 is 
|p| = bee 

| (2 Zo)? + (AX2)? 
f AB, is made equal to zero at 350 mc (0.058 wave- 
ength stub) by the use of a suitable fixed inductance, 
he reflection coefficients of Fig. 4 would result for the 
tub of measured By given in Fig. 3. 


|p| = 


(9) 


(10) 


REFLECTION COEFFICIENT | (| 


04 06 08 10 
LENGTH OF STUB (WAVELENGTHS) 
Fig. 4 


a 


J 


Similar measurements made with coaxial-line tuning 
lements and a small loop antenna showed even more 
ronounced effects of losses on the tunable band be- 
ause of the correspondingly smaller radiation resist- 
nces. Such loops have one advantage over dipoles, 
owever, in that the required capacitive tuning ele- 
nents are superior in many instances to dipole tuning 
nductors with regard to losses. 

The tuning efficiency for the circuit of Fig. 1 (with 
oss components added) is obtained by writing the ex- 
ressions for power input and power delivered to the 
ntenna in terms of the circuit voltages and currents. 


‘he result is 
yer cent En 100 


E Cesar 
i 14242 [G+ G+ Bot By! 


a a 


Es ake 
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which for high-Q tuning elements (G;«B,, G.<Bo, and 
G,<G,) and with condition (b) satisfied, reduces simply 
to 


<8 100 
per cent Hy, = 12 
page (12) 
1+ —+ — 
Ga a 
or, for the dual circuit of Fig. 2, 
100 
percent Ly, = (13) 
RR, R, 
1+—4+— 
APPLICATION 


The application of these single-control tuning circuits 
is limited by losses in the tuning elements and in the an- 
tenna. Comparison of those losses with the antenna 
radiation conductance or resistance will indicate the 
efficiency [according to (11)] and the tunable band 
[from (7) and (10)] to be expected. The circuit of Fig. 1 
where high-Q lumped capacitors can be used would 
seem to be most promising. The tunable cavity antenna 
mentioned previously embodies this feature. A tunable 
band of two to one with a reflection coefficient less than 
0.15 has been obtained with this antenna for an aperture 
between one eighth wavelength and one fourth wave- 
length square.® 

The problem of obtaining suitable inductances for use 
in the dual circuit of Fig. 2 is more formidable. The ob- 
vious procedure in this case is (a) obtaining the maxi- 
mum Q inductance for Xj, (b) use of dipole with as large 
R,/Xa ratio as possible, and (c) obtaining an optimum 
tunable band by adjusting the dipole length to achieve 
a compromise between condition (b) and the length for 
a minimum R,. 

The percentage bandwidth as defined above is a prop- 
erty of the antenna alone. It can be maximized by ad- 
justment of the antenna for maximum R, and minimum 
XG. 

Analysis of the circuits was based on the character- 
istics of lumped elements. The desired frequency de- 
pendence of the fixed tuned element, however, also re- 
sults approximately for the case of a short length of co- 
axial line. Such an element may be useful when a fixed 
inductor or capacitor at high frequencies is desired. Co- 
axial lines may also be used to obtain the variable re- 
actances. However, loss resistances and discontinuity 
susceptances must be taken into consideration with care 
in the selection of an element for a particular frequency. 
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Design of Line-Source Antennas for Narrow 
Beamwidth and Low Side Lobes* 


TA DaeTAY LOR | 


Summary—It is well known that the phenomenon of radiation 
from line-source antennas is very similar to that of the diffraction of 
light from narrow apertures. Unlike the optical situation, however, 
antenna design technique permits the use of other-than-uniform dis- 
tributions of field across the antenna aperture. Line source synthesis 
is the science of choosing this distribution function to give a radiation 
pattern with prescribed properties such as, for example, narrow 
angular width of the main lobe and low side lobes. In the present 
article the mathematical relationships involved in the radiation 
calculation are studied from the point of view of function theory. 
Some conclusions are drawn which outline the major aspects of 
synthesis technique very clearly. In particular, the problem of con- 
structing aline source with an optimum compromise between beam- 
width and side-lobe level (analogous to the Dolph-Tchebycheff 
problem! in linear array theory) is considered. The ideal pattern is 
cos 71/u2— A?, where u=(2a/)) cos 0, a is the half-length of the 
source, and cosh 7A is the side-lobe ratio. Because of theoretical limi- 
tations, this pattern cannot be obtained from a physically realizable 
antenna; nevertheless its ideal characteristics can be approached 
arbitrarily closely. The procedure for doing this is given in detail. 


INTRODUCTION 


HE term line source may be applied to any an- 

tenna which embodies a long, narrow, and straight 

geometry and which depends upon variations in 
field or current strength with respect to the longitudinal 
co-ordinate as the basis for its directivity. It is implied 
that these currents or fields are continuous functions of 
this longitudinal co-ordinate, hence arrays of discrete 
elements will not be considered as such. It is well known, 
however, that arrays of a large number of elements are 
most conveniently analyzed by regarding them as con- 
tinuous distributions, therefore the remarks made here 
will apply—in this sense—to such arrays. 

The pattern of a line source contains an element factor 
and a space factor. The element factor depends upon 
the type and orientation of fields or currents which 
make up a typical segment of the source, and usually 
has little directivity, at least in the direction in which 
it would augment the directivity of the space factor. 
The space factor is the highly directive factor; it is 
dependent for its shape upon the relative variation of 
field or current strength along the source and is the ob- 
ject of the antenna designer’s art. In shaping the space 
factor, the designer may or may not wish to take cog- 
nizance of the element factor. In this article it will be 
assumed that he does not and that the objective is to 


* Original manuscript received by the PGAP, April 26, 1954; 
revised manuscript received October 11, 1954. The work reported 
in this article was performed at the Hughes Aircraft Company and 
has been described in Hughes Aircraft Company Tech. Memo. No. 
316. The work was sponsored by the Air Force Cambridge Research 
Center, Cambridge, Mass., under Contract AF 19(604)-262. 

} California Institute of Technology, Pasadena, Calif.; formerly 
Research Physicist, Hughes Aircraft Co., Culver City, Calif. 

'C. L. Dolph, “A current distribution for broadside arrays which 
optimizes the relationship between beamwidth and side-lobe level,” 
Proc. I.R.E., vol. 34, pp. 335-348; June, 1946. 


design a space factor with the narrowest beamwidth 

and lowest side lobes consistent with both the available 

aperture and the necessity for avoiding a super-gain line 
source. 

Let the line source, which has a total length 2a, have 

a distribution function g(x). Then the space factor is 

given by the familiar Fourier-transform relationship 

a 
S(v) = ef g(xje*"dx, (1) 
where vy=cos @ as illustrated in Fig. 1 and k is the free 


space wave number, 27/d. It will be advantageous to 
introduce new variables as follows: 


2S = (2) 
a 

= ta (3) 
u = x 


The quantity \/2a will be known as a standard beam-. 
width? and will prove useful later. With these substitu-: 
tions, (1) becomes 


20.0 
sw) = = if o(peudp. (4) 


The factor 2a/d will be dropped and the integral de-. 
noted F(u): 


Fu) = [s(perap. 


aes 


(5) 
Eq. (5) expresses the basic relationship which will be: 
studied here. . 


LINE SOURCE 


Fig. 1 


PROPERTIES OF THE DISTRIBUTION FUNCTION } 


It now becomes necessary to specify the hypothese 
which the function g(p) must satisfy in order that th 
results of this analysis will be applicable. An exhaustiv. 


* The beamwidth between half-power points of a uniformly il 
luminated line source is 0.8859 of a standard beamwidth. 
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snerality is not attempted here; nevertheless, it is felt 
iat a high degree of relevancy to the physical situation 
achieved. 

It is necessary, in antenna work, that g(p) be single- 
alued, uniformly bounded, and at least piece-wise con- 
nuous in —7 <p<7. However, anyone conversant with 
ich work will grant that there exists a large class of prac- 
cally important cases in which not only g(p), but all its 
erivatives as well, are continuous over the full in- 
erval —7<p<m. This article is specifically concerned 
ith functions of the latter type and all general state- 
rents made henceforth are to be understood as apply- 
ig only to distributions which obey this hypothesis. 
he restriction, however, is not as stringent as might be 
upposed. All analytically expressible functions which 
re non-singular in the aforementioned interval, such 
s the uniform function cos p/2, Jo(p) (and an infinitude 
f others which have no specific names), belong to this 
lass. In fact, any function which is continuous in the 
iterval under discussion can be approximated arbi- 
rarily exactly by a member of this class. This article 
rill also, incidentally, be concerned only with even dis- 
ribution functions, since these are sufficient for the 
urpose at hand, namely the optimization of the com- 
romise between beamwidth and side-lobe level. 

The degree of abruptness with which the source dis- 
ribution begins and ends at p= +7 will be a matter of 


onsiderable importance. Let the symbol “~” be read 
is asymptotic to.” Then suppose that 
eee ey a) as p> 1 (6) 


g(p) ~ K(r — p)* as pot) 


‘he quantity K is a non-zero constant in these equa- 
ions. The real parameter a, if an integer, is equal to the 
aultiplicity of the zeros in g(p) at p=+7. Thus if 
-=0, g(p) assumes the fixed non-zero value K at these 
wo points and the distribution is said to have a “pedes- 
al.” If a=1, the distribution falls to zero linearly at its 
nds as in cosine distributions. If ~=2, it falls to zero 
uadratically as in cosine-squared distributions. All 
alues of a greater than —1, both integral and non- 
itegral, will be considered. The not-uniformly-bounded 
istributions corresponding to a in the range —1<a<0 
re included in order to achieve a degree of completeness 
ecessary for later arguments, even though such values 
re not regarded as being physically realizable. If g(p) 
. known analytically, a may be found by taking the 
mit as p approaches —7 of (p+7) g’(p)/g(p), or the 
mit as p approaches 7 of (pb—7) g’ (p)/g(p). Evidently, 
‘is possible to write g(p) as the product 


g(p) = h(p)(m? — p*)*, (7) 


here h(p) is an even function which does not vanish at 
= +7 and which fulfills the other requirements apply- 
ig to g(p). 

Let us now reconsider (5) from the point of view of 
inction theory. Let the variables p and u be imbedded 
the complex domains & and z respectively, such that 
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E=p+ig and z=u-+iv. The functions g(p) and F(u) 
are then profiles on the axes of reals of the functions 
g(&) and F(z). The visible range of the z-plane® is the 
real segment —2a/ASuS2a/X, and a unit distance in 
the z-plane corresponds to one standard beamwidth 
(see Fig. 2). The integral in (5) becomes the limit at 


VISIBLE RANGE 


jet: 


Z PLANE 


Fig. 2 
7 tends to zero of the integral along C of Fig. 3, begin- 
ning at —a-+7 and following the axis of reals to t—7 


in the & plane. The limiting process is necessary because 
of the formal consideration of negative values of a. 


iq 


Ce intone mane) . 


€ PLANE 
Fig. 3 
Eq. (5) becomes 
F(z) = lim | g(é)e*dé. (8) 
™0 (6) 


This g(€) is, of course, equal to the product A(é) (7? — 
&)«, If we cut the — plane from — © to —z and\from 
+77 to +, then (7?—£&?)¢ is regular in this cut plane 
and should be regarded as being on its principal branch. 
Since it is desired that the asymptotic behaviour of 


3 An observer who travels around the line source to examine its 
radiation pattern will cover all real angles in the range 0 <0 <7. This 
coverage is translated into the point set —2a/ASuS2a/d, v=0, in 
the z-plane by the relationship “= (2a cos @)/ as shown graphically 
in Fig. 2. Hence the term visible range is applied to this point set. 
The profile of F(z) on the visible range is the only part of F(z) which 
ormally enters into the formation of the radiation patterr- 
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g(p) as p approaches +7 be as indicated in (6), 1e., 
be embodied entirely in the factor (7?—£)¢, it is essen- 
tial that h(£) be regular at least in neighborhoods of 
—r and +7 and be different from zero at these points. 
Moreover, if g(p) and all its derivatives are to be con- 
tinuous in —7t<p<z, h(&) must be regular in a domain 
which contains this interval in its interior. Thus it is 
sufficient if h(£) is different from zero at £=+7a and 
regular in the race-track shaped region R of Fig. 3, 
where 6 is the arbitrarily small but non-vanishing radius 
of the ends of the region. 

The hypotheses of this analysis may therefore be 
summarized as follows: 


1. g(p) is the profile on the real interval —7t@<p<m 
of the principal branch of the function A(&) (7? — 
ENS 

. h(€) is regular in the region R of Fig. 3. 

. h(&) does not vanish at E= +7. 

. h(&) is even. 

. ais real and must be greater than or equal to zero 
for physically realizable distributions; however, 
values of a between —1 and zero will be considered 
for completeness of analysis. 


On Bm W bd 


PATTERN THEOREMS 


In this section, some basic theorems concerning the 
properties of the function F(z) will be proved. 


Theorem I: F(z) 1s an even function of z. 


This is a very simple theorem. Because of the even- 
ness of A(&) and of the factor (7?—£&)¢, the transforma- 
tion may be written as follows: 


F(z) = lim eee 2h(p) (x? — p)* cos padp, (9) 
™0 


0 


since cos p(—2) =cos pz, F(—z) = F(z), and F(z) is even. 


Theorem II: F(z) is an entire function of z. 


Eq. (9) will be adopted as the point of departure in 
the proof of this theorem. It will be necessary and suffi- 
cient to show that: 


1. The integrand is itself entire for every p in the 
interval OS pS7—-r. 

2. The integrand is continuous with respect to p for 
every p in the interval 0<pSXm-—7 and for every z 
in the entire finite plane. 

3. The limit of (9) is approached uniformly for g in 
G, where G is any closed bounded subset of the 
entire finite plane. 


Conditions 1 and 2 are evidently satisfied. To demon- 
strate the fulfillment of condition 3, it must be shown 
that the difference, f(z, 7), between the integral in (9) 
and its limit is uniformly bounded for zg in G and that 
this bound tends to zero as 7 tends to zero. Evidently, 


f,2) = [2H pV + pam — p)* cos ped. (10) 
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A function is always uniformly bounded on a closed 
bounded set within its domain of regularity. Thus, for 
p in the closed interval r—7 S$pSz and z in G, | 2h(p) 
(r+p)¢| is bounded by some constant A, and since 
cos pz is entire with respect to z for every value of p in’ 
the stated interval, | cos pz| is bounded by another 
constant B. Therefore 


[fe |S 4B f (—pyrap (11) 
A Brito 
OD pete nes 3 12 
lf, 7)| S ey (12) 
Since a>—1, 
lim | f(z, 7) | = 0 (13) 


7™0 


for every z in G, hence condition 3 is satisfied and F(z) 
is entire. For further background in connection with 
this proof the reader should consult Copson.* 


iq 


€ PLANE 
Fig. 4 


Theorem III: as z tends to infinity, F(z) has the following 


asymptotic forms 
Gree) 
cos r{| 2 — 
2 


1(z)1+e 


for Rez >0O (14) 


(oie) 
COS = — 
2 


1( = 2) l+a 


for Re 2:<.0,, (45) 


F(z) ~ h(r)(2r)**9T(1 + @) 


F(z) ~ h(r) (20) +T(1 + @) 


etlzl 


Bs) Ma) OS ae 


for Rez =.0. (16) 

Proof of (14) depends upon deforming the contour C 
into the contour C, C, C; and subsequently allowing the 
radius of the small arcs to shrink to zero (see Fig. 4). 


*E. T. Copson, “Theory of Functions of a Complex Variable,” 
Oxford University Press, pp. 107 ff., 116; 1935. 
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It is easily shown that even for —1<a<0 the con- 
tributions along the small arcs vanish. We begin by 
saying: 


= re" 


(17) 


Tv Tv 

Chase OUST CONStAM tre 
y D 

Ti 


The angle ¢ determines the direction taken by the two 
parallel legs of the contour. The length of these legs is 
b, 7<b<6. Let J; be the integral along C,. Then 


é = — TT — pe<7/2—) 
dé = dpe'!/2—9) | (18) 
b 
ic = lim h(&) (4 she £)“preia(n/2 d) @ ine prei(m/2 ®dp (19) 
70 T 
b 
Men ett @)(x/2-9) Jim f h(&)(a — £)%p%e-*dp. — (20) 
tT 0 T 


Make a change of variable. Let p=¢/r: 


em (2—(C+a) /2) br te?(t/ 2-9) 
= ———————. lim h\) — 7 + ———— 
: gita 7T0 TT. | if 
tei(7/2—4) a 
| 2 — | bedi: (21) 
A 


The only significant contribution to this integral takes 
place at some finite ¢, specifically at #=0 if aSOor at 
t=a if a>0. Hence two of the factors of the integrand 
can be removed as ¢ tends to infinity; these become 


simply i(—7) (27)*. Thus: 
br 
ff bi Cmte (22) 


e—ix(2—G+a)/2) h( — 7) (27) ; 
h= lim 
gita 70,7” 

Since 7 goes to zero independently of 7, we may substi- 
tute zero for the lower limit. Then as 7 tends to infinity, 
the upper limit tends to infinity and the integral tends 
to (ita). Hence, 


e7 tt (z—(1+a)/2) 


TI, = h(—7)(20)**T(1 + a) ante) (23) 
Similarly, 
eit (e—(1t+@)/2) 
Iz = h(mr)(2r)'*eT(1 + a) On (z) ite (24) 
On C; we have 
ay i(m/2—¢) 
& = p+ be | (28) 
dt = dp 
l= f Wee Byrememap, (26) 


Here the function h(£) (7?—£)* is uniformly bounded 
on C2 by some constant B. Then 


| Is| S 2rBerre™ |¢ |. (27) 
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But the magnitude of one or the other of the two inte- 
grals J, or J; (depending upon the sign of #) tends to 
TGP Ne" el orelarge 7.) | heretore | Z| will, for 
large r, always be negligible with respect to | 7, or 
| 73| or, if ¢6=0, to both | 7 and | Zs]. The total inte- 
gral is thus the sum of (23) and (24). Since h(—7) = 
h(a) this sum is simply (14). Exp. (15) may be derived by 
deforming the contour as in Fig. 5, or by appealing to 


iq 


€ PLANE 


Fig. 5 


Theorem I. Exp. (16) is simply the limiting form of both 
(14) and (15), and can also be derived by the methods 
used in obtaining (14) and working on the undeformed 
contour C. It is interesting to notice that when a is a 
non-negative integer there are no singularities at = +7. 


TABLE I 
a asymptotic form of F(z) 
0 Mn) 
12 
; i(m)(2n)? COS 72 
12? 
2 h(w)(2n)(2) 
TZ 
3 h()(2m)*(6) — 
1 


The contour of Fig. 4 can be swung around until it con- 
forms to that of Fig. 5 without making a separate an- 
alysis. The net result, that one asymptotic expression 
suffices for all g, | z| sufficiently large, is shown in Table 
I above. As a further consequence of Theorem ITI, it is 
seen that the asymptotic behaviour of F(z) depends only 
upon the manner in which g(p) acts at the ends of the line 
source. If there is a pedestal, the asymptotic form must 
be sin 12/72. If g(p) falls to zero linearly at the ends of 
the line source the asymptotic form must be cos 12/72?; 
if quadratically, it must be sin wz/mz*, etc. The decay 
rate of the remote side lobes in F(z) is thus irrevocably 
fixed by conditions at the ends of the line source. 
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Theorem IV: Let the zeros of F(z) be grouped in pairs: 
+21, +22, etc., and numbered in such a way that the modult 
of the members of successive pairs do not decrease. Then 
the members of the nth zero pair tend to the positions + 
(n+a/2) as n tends to infinity. 


La 


PLANE 
Fig. 6 


In Fig. 6, circles of very small fixed radius € have been 
drawn around the points +[T+(a/2)+k]| where T 
and k are integers and k=1, 2, 3, -- +. The large circle 
H has a radius of T+(a/2)+(1/2). Now if T is made 
large enough but finite, this system of circles will cover 
all the zeros of F(z), for if the circle H is sufficiently 
large, F(z) will approach its asymptotic form arbi- 
trarily exactly in the region outside of H and the zeros 
of the asymptotic form of F(z) do lie at the points 
+ [T+ (a/2)+k]. The question which must be settled 
is that of the number of zeros within H. These zeros 
constitute a bounded set and their number cannot be 
infinite, for if it were they would have a limit point and 
F(z) would have an essential singularity within H, a 
condition which contradicts the entirety of F(z). Let 
the number of these zeros be 2M, M an integer. Then 
F(z) can be approximated arbitrarily exactly by the 
product of a polynomial of degree 2M and a canonical 
product on zeros outside of H located at the points 
+ [T+(a/2) +k] where & runs from one to infinity.® 

The only defect in this approximation is embodied 
in the small distances between the actual zeros of F(z) 
and the points +[T+(a/2) +k] but these distances can 
be made smaller than any fixed ¢ by taking 7 large 


5 The technique of synthesizing the space factor of a line source 
by setting up the product of a polynomial and an entire transcen- 
dental function where the latter has the same zero arrangement as 
that just discussed here has been mentioned earlier by the author in 
less general terms in an unpublished paper, “An antenna pattern 
synthesis method with discussion of energy storage considerations,” 
presented at the U.R.S.I. Washington Meeting; April 1951. Here, the 
quantity T+(a/2) is equal to the N/2 in that paper; only cases cor- 
responding to non-negative integer values of a were treated. 
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enough. The approximation, then, is 


F(z) = Paw(s)[]]1— e (28) 


a ae aC 
ms (7 Si arb t) 
2 
This canonical product is expressible in I’ functions, for 
a 2 
(r+ 544) 


(a4 a 
r(r+S4i14s)r(r+ 541-2) 


a) 32 i 
=p E ~ z mar 
il (z +—+ z) 
Z. 
Hence from (28) the asymptotic form of F(z), as | z| 
becomes much greater than the radius of the circle H, is 


Kz?! 


(29) 


F(2)~ (30) 


a a 
r(r+S414s)r(7+ 241-2) 


where K is a constant including both the coefficient of 
2™ in the polynomial and [I(T +a/2+1) ]?. Now 


1 sin 7A 
= . (31) 
I(A)P(1 — A) T 
Let A = [z—T—(a/2)]. Then 
sine (< _ <) 
2 
F(z) ~ (—1)?K224%@ ————___——_. 
r(s —T- <) 
2 
(32) 


r(s+7+5+1) 


By appealing to the Stirling approximation for I'(u), 
valid where |u| is large and where a wedge-shaped 
region covering the negative axis of reals is excluded, 
it can be shown that 


Qa 
r(:-7-<) \ 
2 1 


eee os (33) 


2 T+1+a 


r(s+r+241) 


in a region which is even more extensive than that de- 
fined by Re z>0. Then certainly for Re z>0 


1 
(—1)7Kz™™ cos r(s a =) 


Egy (34) 


oeTygita 


But by Theorem III, F(z) is asymptotic to a constant 
X cos t[z—(1+a)/2]/rz't# in the region defined by 
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Re z>0. Hence we must conclude that M=T and that 
there are 2T zeros within the circle H. Therefore the zero 
within the € neighborhood of T7+(a/2)+1 must have 
been a member of the T+1st pair of zeros, that within 
an € neighborhood of T+(a/2)+2 must have been a 
member of the T+2nd pair, etc. Therefore zg, tends to 
n+(a/2) as n tends to infinity, and the theorem is true. 


Basic DESIGN CONSIDERATIONS 


We are now in a position to state some definite con- 
clusions regarding the design of line sources for highly 
directive patterns. We have seen that the pattern of a 
finite line source, the even distribution function of which 
decays with the exponent a at the ends of the source, 
must have the following essential characteristics: 


1. F(z) must be an even entire function of z, z=u-+iv. 

2. The remote side lobes of F(z) must decay as 
| «| —Gt@), 

3. The remote zeros of F(z) must lie on the axis of 
reals and the members of the mth zero pair must 
approach the positions + [7+(a/2) | as tends to 
infinity. The remote zeros are separated by unit 
intervals on the real axis. 


Since, from the function-theoretic point of view, two en- 
tire functions with the same asymptotic behavior are 
_ but trivially different, we may expect any entire func- 
tion possessing the characteristics noted to be the space 
factor of some hypothetical line source. Space factors 
will differ from one another, then, only in two aspects: 


1. The value of a. 
2. The placement of the central zeros; i.e., the zeros 
near the origin in the z-plane. 


As with the space factors of arrays of discrete ele- 
_ ments, zeros are precious items. In highly directive pat- 
terns they are used to “hold down” the average level of 
F(z) in the side-lobe region and their absence over even 
a relatively short range permits a large lobe to appear. 
As an aid in visualizing this, the harmonic function log 
| F(z) | = U(u, v) may be regarded as describing an elastic 
surface which is supported at infinity and weighted 
down by small heavy weights at each of the zeros. 

Since we are primarily concerned with the profile of 
F(z) along the axis of reals, it is important to realize that 
the average value of this profile in the side-lobe region 
can best be minimized if all the zeros are on the axis of 
reals. A zero just off this axis produces a dip but not a 
null in the profile. Furthermore, it is a matter of expe- 
- rience® that if, in a given portion of the side-lobe region, 


6 As an example, consider the entire function fi(z) =cos 72, which 
has equally spaced zeros and uniform lobes. Move the two zeros 
which straddle the origin together to form a double zero at the 
origin. The function becomes 

2? COS 72 


fo(z) = eon ts ae x 


The remote behavior is unchanged but the two central lobes of fo(z) 
have risen to 1.38 compared with 1.00 for fi(z). 
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a fixed number of zeros are available, the average side- 
lobe level will be lower if these zeros are spread out 
approximately uniformly rather than grouped in twos 
or threes or some other multiplicity. Therefore the first 
design conclusion is that for highly directive patterns, 
the central zeros of F(z) should be simple zeros on the 
axis of reals. The remote zeros must follow this law 
anyway, hence, for highly directive patterns all the zeros 
of F(z) should be simple zeros on the axis of reals. 

Let us next consider the selection of a value of a. 
First of all, to insure that g(p) will be uniformly bounded 
a must be greater than or equal to zero. But as a is in- 
creased the zeros of F(z) abdicate from the central 
region of the z-plane, since the uth zero pair tends to the 
position +[n+(a/2)]. Increasing a by one unit is 
tantamount to throwing away one zero. It is true that 
high values of a mean rapid decay rates of the remote 
side lobes, but since any side-lobe structure within the 
central region can be obtained anyway, and since the 
central region can be extended to cover the visible range, 
—(2a/d) Su S(2a/d), a rapid decay rate of the remote 
side lobes will not be a normal desideratum in line- 
source design. It follows, then, that a should have the 
smallest possible value, zero, and that the best distribution 
function for generating a highly directive pattern will have 
a pedestal. 

If @ is to be zero, the remote zeros must be asymptot- 
ic to the integers; the whole design technique becomes 
a question of positioning the central zeros along the 
axis of reals in such a way that the best beamwidth is 
obtained compatible with the level of the central side 
lobes and in such a way that the members of the nth 
zero pair smoothly approach the positions +n. 


AN ENTIRE FUNCTION WITH UNIFORM SIDE LOBES 


Before proceeding with the technique of practical 
line-source design, it will be essential to formulate the 
ideal space factor, i.e., an entire function with uniform 
side lobes. Later it will become apparent that a practical 
line source can have a space factor whose characteristics 
approach the characteristics of this ideal space factor 
arbitrarily closely and that the ideal cannot be imitated 
exactly without resorting to super-gain. 

G. J. van der Maas’ has shown that as the number of 
elements of a Dolph-Tchebycheff array is indefinitely 
increased, the space factor approaches a function which 
(expressed in the present notation) is 


Fo(u, A) = cos r-V/u? — A?. (35) 


Here A is an adjustable real parameter having the prop- 
erty that cosh 7A is the side-lobe ratio. This function 
(35) is evidently the profile on the axis of reals of the 
entire function 


F(z, A) = cos r+/2? — A’, (36) 


7G. J. van der Maas, “A simplified calculation for Dolph- 
Tchebycheff arrays,” Jour. Appl. Phys., vol. 25, pp. 121-124; Janu- 
ary, 1954. 
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which can be shown to have the zeros 


Zn = + aA/A2+ (wn — 4)?; w= 1,2,3--- @. (Sin) 
This entire function is clearly the ideal space factor. 
Its side-lobe region (w?> A?) is populated by an infinite 
number of equal lobes of unit amplitude; its main lobe 
appears in the region u?<A? and has adjustable ampli- 
tude, as was pointed out. Nevertheless this space factor 
is unrealizable beecause the remote side lobes do not 
decay and the members of the mth zero pair are asymp- 
totic to +[n—(1/2)], both of which conditions corre- 
spond to an @ value of —1, a value which contradicts 
the uniform boundedness of g(p). However, any good 
practical space factor must have characteristics resem- 
bling those of Fo(z, 4) in the central region of the z- 
plane, but in the remote regions its side lobes must de- 
cay as 1/z and its zeros must tend to the integers, not 
the integers minus one-half. The problem is becoming 
clearer: in the central region the zero positions must be 
related to those of Fo(z, 4) but in the remote regions, 
the mth zero must tend to +n. A solution to this prob- 
lem is proposed below. 


SYNTHESIS OF A PRACTICAL LINE SOURCE 
APPROACHING ARBITRARILY CLOSELY 
TO THE IDEAL 


Let the practical space factor resemble the following 
function in the central region: 


Z 2 
F(z, A) =<cos r4/(=) — A’, (38) 
oO 
the zeros of which are 
Zn = taoV/A2+ (nm — 3). (39) 


The constant o is a number somewhat greater than 
unity. One can visualize (38) as an ideal space factor 
slightly dilated horizontally. Now let us choose o such 
that eventually one of the zeros of Fi(z, A) becomes 
equal to the corresponding integer, designated #. Then 


f= bn >= £oV/A* +" — 3)%, (40) 


whereupon @ is defined by 


n 


"VAP (HE 


(41) 


We now propose to construct an entire function with 
the following zeros: 


Zn = + V/A? + (nm — 3)? 


ey = ae 


The canonical product is 


n—1 g? 
Fg, Ayn) =C {1 - \ 
I o?| A? + (n- 1)2] 
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(a) 32 
II (1 - *). (43) 
n=n n® 

Let the second product be denoted Qaa—»(z) to agree 
with the notation of the unpublished paper referred to 
earlier. This product has some convenient equivalent 


expressions: 
Qs @—1) (2) = I (1 7 =) (44) 
ren) as 
(av + 2) (% — 2) 
i Coaeses Will (46) 
(7 —1+ 2)!” —1-—=2)! 
a = TZ = 47) 
rz [J (1 = =) 
n=1 n” 


The entire function with the zeros of (42) is 


n=1 


n—1 32 
F(z, A; A) = CQ2a-1(2) I] Sosnpoweaa ° (48) 


Here C is an arbitrary constant; if it is taken equal to 
cosh 7A then F,(0, A) =F(o, A, 7). Let us suppose that 
this has been done, then 

lim F(z, A, m) = Fo(2, A); 


nave 


(49) 


for as # increases, more and more of the zeros of F(z, A, 
n) follow the pattern of (39) and simultaneously o—1, 
hence (49) is as stated. Thus F(z, A, v) defines a two- 
parameter family of entire functions whose members ap- 
proach the tdeal space factor arbitrarily closely as n is in- 
creased. The value of #, which is finite in any practical 
case, is of considerable significance. In relationship to 
the pattern, the points + 7 divide the region of uniform 
side lobes (| | <*n) from the region of decaying side 
lobes (| u| >) where here, as before, z=u-+iv. As fi is 
increased, the region of uniform side lobes extends far- 
ther out from the main beam. In the region of uniform 
side lobes the lobes actually decrease slightly as | 2e| 
increases and their initial level is very slightly lower 
(that is, better) than the design level, 7. (Note that 
mA =arc cosh 7.) The chief disparity between the prac- 
tical pattern and the ideal is in the beamwidth, which is 
greater than the ideal by a factor almost exactly equal 
to o. However, # does not have to be very large to make 
a only a few per cent greater than unity.’ 

To summarize, the members of the pattern family 
F(z, A, x) have two independent characteristics: 


1. The design side-lobe ratio, 7. 
2. The boundary of the region of uniform side lobes, 
n, an integer. 


* For example, if the design side-lobe level is to be 25 db, » =17.78; 
A =1.137. Then if m=5, ¢=1.077 and the beamwidth of the practical 
pattern will be only 7.7 per cent greater than that of Fo(z, 1.137) 
the ideal pattern. : 
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Q3 
TABLE II 
(1) (2) (3) (4) (5) 
ae 5 18085 
esign side- n r S e 
eee (Side tobe . Be Values of the parameter o 
(in db) voltage ratio) | (in deg.) a) @=e n=4 n=5 n=6 a n=8 
: 1.00000 28.65 0.00000 il Sees) 1.20000 1.14286 PATI 1.09091 1.07692 1.06667 

i 1 peas 34.49 0.14067 1.29351 1.18672 1.13635 MOU 1.08838 1.07514 1.06534 
ty 3.16228 40.33 0.33504 1.24393 1.16908 TSA: 1.10203 1.08492 1.07268 1.06350 
- 5.62341 45.93 0.58950 1.18689 il A 1.11631 1.09528 1.08043 1.06949 1.06112 
20 10.0000 Sul silky 0.90777 1.12549 ieied3ss iL 1O273 1.08701 1.07490 1.06554 1.05816 
a ce 56.04 1.29177 _ 1.09241 1.08698 1.07728 1.06834 1.06083 1.05463 
a 31.6228 CORSON seem 74220 — — 1.06934 1.06619 1.06079 1.05538 1.05052 
B.S 56.2341 64.78 2.25976 — — —- 1.05386 1.05231 1.04923 1.04587 
40 100 .0000 68.76 2.84428 — _- — 1.04298 1.04241 1.04068 


Once these two parameters have been chosen, the pat- 
tern, the distribution function, and all other relevant 
data may be calculated. In selecting , it is essential to 
avoid values that are too small; the change from F;(z, A) 
to F(z, A, v), in which the trans-7 zeros migrate to the 
integers, must be such as never to increase the spacing 
between any of these zeros. In order that this condition 
be fulfilled, % must be chosen such that a unit increase 
in % does not increase o. Practically, this means that for 
a space factor with a design side-lobe ratio of 25 db, % 
must be at least 3, and that for a design side-lobe ra- 
tio of 40 db, it must be at least 6. Apart from this lim- 
itation one has considerable liberty in the choice of 
n; large values make o more nearly equal to unity, there- 
by sharpening the beam; however, the additional benefit 
obtained in this way soon becomes negligible. If 7 is 
increased sufficiently beyond the endpoint of the visible 
range, 2a/X, all the visible side lobes tend to become uni- 
form. Such large values of 7% have the effect of super- 
gaining the line source, a phenomenon which will be 
considered in detail in the next section. The choice of 
# may also be influenced by considerations of gain. 
' The beamwidth of the ideal space factor (in standard 
beamwidths) is a function of the side-lobe ratio only. 
It may be derived® from (35) and is given rigorously as 
follows: 


2 4/ ( h )2 ( h n ) 
aaa a arc.COsSh ——] © 
‘ec arc cosh 7 iW 


A curve showing the relationship (51) is presented in 

Fig. 7. The beamwidth of an F(z, A, #) space factor is 

given by the following extremely good approximation: 
AB 

Ja + (i — 


(50) 


B = oBo = (51) 


Table II (above) shows relationships among various 
quantities discussed here. Column (1) of the table con- 


-~9R. J. Stegen, “Excitation coefficients and beamwidths of 
Tchebycheff arrays,” Proc. I.R.E., vol. 41, pp. 1671-74; November, 
1953. 


Stegen has derived essentially the same formula by considering 
Dolph-Tchebycheff arrays of large numbers of elements. 


: 


tains the design side-lobe ratio in decibels—a quantity 
which is usually chosen first in any synthesis proce- 
dure. Columin (2) is self-explanatory. Column (3) gives 
Go expressed in degrees. This must be multiplied by 
d/2a to get the ideal beamwidth for a given line source. 
Column (4) gives the parameter A*®. Column (5) gives 
the parameter o for several possible choices of 7%. As 
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Fig. 7—Beam width vs. side-lobe ratio for ideal space factor. 


an example, suppose a line source with a 30-db side- 
lobe level is required and that the aperture-to-wave- 
length ratio is 50, that is (2a/A) =50. Then the ideal 
beamwidth is 60.55/50 or 1.211 deg. If it is decided to 
use an # of eight, ¢=1.05052, and the practical beam- 
width, that is, the beamwidth of the function F(z, A, n), 
is 1.272 deg. The zeros of the pattern and the pattern 
itself may be readily calculated using o and A’. The 
places in the table where there is no entry for ¢ corre- 
spond to choices of # which are too small, as has been 
explained earlier. 
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GAIN AND SUPER-GAIN IN LINE SOURCE ANTENNAS 


In this article the term gain of a line source means the 
ratio of power flux radiated in the direction of maximum 
intensity to the value of this flux averaged over all di- 
rections. If 6, @ are spherical co-ordinates in terms of 
which the radiation field is described and if P(@, 6) is 
the element factor (with dimensions of power-per-unit 
solid angle), the total power flux pattern is simply 
P(O6, >) | S' (cos 8) | 2, With 6,, 61 representing the direction 
of maximum intensity, the gain is given by 


re P(61, $1) | S(cos 61) 2 - (52) 


1 2a T 
— f f P(6, ¢) | S(cos 0) |? sin 6d6dp 
dr Jo 0 


It has been mentioned that P(6, ¢) usually has little 
directivity, especially in the 6 co-ordinate. If S(cos @) is 
strongly directive and not endfire (6; not near 0 or 7), 
then a very good approximation for the gain is obtained 
by assuming that S has the properties of an impulse 
function and that P(61, ¢) can therefore be factored out 
of the @ integral 


P(A, 1) 


Qr 
if P(6,, ¢)d6 


G= 
1 
ar 0 


| S(cos 6) [2 


ae (53) 
—{ | S(cos 6) |? sin 6d0 
2/0 


The second factor in the above will be called Gp and is 
the gain which the entire antenna would have if the 
element factor were isotropic; from the point of view of 
line-source design, Gp is the quantity which one would 
naturally discuss when comparing various line sources 
on the basis of gain. In the present case the beam is 
presumed to be broadside and cos 0,=0. Replacing cos 
6 by v, Go becomes 


S(0) |? 2 
eoeene Oy mee On Cer sree. 
1 1 : 2a/d : 
mie | S(v) |2dv J Fo [raw 


The quantity \Go/4a, which will be called the specific 
gain, is of considerable importance in what is to follow. 
It is, of course, given by 


| FO) |? 
2a/r 


| F(u) |2du 


—2a/r 


AG) /4a = (55) 


Now consider a line source with a distribution g(p) of 
fixed functional form and let the ratio of physical length 
to wavelength be increased. This produces no change in 
F(u), which is related to g(p) through (5), however, it 
obviously changes the specific gain via the limits of the 
integral in (55). For very small values of 2a/\, the 
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specific gain changes rapidly with respect to this quan- 
tity; when, however, 2a/\ passes into the region where 
the side lobes of F(w) have become low and are getting 
lower by virtue of their asymptotic decay properties, 
the specific gain becomes practically constant and un- 
dergoes but little further change if 2a/) increases all the 
way to infinity. As an example, the specific gain of a 
uniform distribution is plotted with respect to 2a/X in 
Fig. 8. Evidently, the limiting value of the specific gain 
is an important line-source parameter, especially since 
it depends only upon the functional form of g(p). This 
is shown by the following expression, in which the Par- 
seval formula’ for Fourier integrals has been employed: 


| f (Pap 


an f "| (p) [tap 


2 


| 7 P 


(AGo/40) 0 = (56) 


foo} 


{h F(u) |2du 


% 


dG,/4a for uniform distribution or 


If g(p) is now varied, it is easy to prove that (AG)/4a),, 
takes on its greatest possible value, namely unity, when 
g(p) is uniform. It is in this sense, and only in this sense, 
that a uniform distribution has a specific gain greater 
than that of any other type of line-source distribution. 

When 2a/X is finite, distributions with arbitrarily 
high specific gains can be synthesized—at least on 
paper. Starting with the space factor of the uniform 
distribution, sin 7u/mu, one can form the space factors 
for such super-gain distributions merely by rearranging 
the existing zeros to give an increased concentration of 
zeros in the visible range. This larger supply of zeros 
can be used to reduce the width of the main lobe while 
holding the relative level of the visible side lobes con- 
stant, thus enhancing the specific gain. At the same time, 
however, (AG)/4a),. must decrease because it had its 


‘0 E. C, Titschmarsh, “Theory of Fourier Integrals,” Oxford Uni- 
versity Press, pp. 50-51; 1948. 


V55 


reatest possible value before the rearrangement of the 
eros took place. It follows that super-gaining must in- 
reaseé the value of the ratio of NGo/4a to (NGo/4a).. with 
2spect to the value it would have if the distribution were 
niform. This ratio, which has been called the super- 


ain ratio,’ is given by 
F(u) |?d 
AGo/4a ia | (w) | # 


o- (AGo/4a) ra {i 


(57) 
| F(w) 2d 


—2a/d 
Vhen this ratio is evaluated for a uniform distribution, 
t will be designated yo; the expression for this is 


w/2 
sin? (27a/d) 
(2ra/h) 


ince (AGo/4a)..=1 for a uniform distribution, yo is 
jumerically equal to the specific gain in that case and 
s therefore also illustrated by Fig. 8. 

It is seen from (57) that y is the ratio of the total area 
o the visible area under | F(u) | *; evidently, then, super- 
‘ining must increase the relative level of the invisible 
obes. Large values of y also imply large ratio between 
nean square value of g(p) and visible area under | F(u) | 2 


1 Tv 
Fe ibe | g(p) |2dp 
as ; (59) 


1 2 2a/r 
(—) i | F(u) |2du 
20 —2a/r 


[his means that, given two line sources with radiating 
ements which are identical and with respect to which 
he approximation of (53) applies, then the mean square 
ralues of g(p) in the two cases will be to each other as 
he two values of y, when both sources are radiating the 
ame total power. This is but one of many ways of 
eaching the conclusion that line sources with an ap- 
yreciable amount of super-gaining are impractical. 

Up to this point, the term super-gain has been used 
o refer to line-source distributions whose specific gains 
re higher than that of a uniform distribution with the 
ame ratio of physical length to wavelength. It has been 
roved that such sources must have associated with 
hem values of y which are greater than Yo. The con- 
rerse is not necessarily true, for high values of y may 
esult from the movement of zeros into the visible range 
or the purpose of giving more character to the pattern, 
s in beam shaping, in ways which would not increase 
he specific gain. The debilities of super-gain sources 
esult, not from the increased gain, but from the in- 
reased value of y. For this reason, it would be prefer- 
ble to extend the definition of a super-gain source to 
iclude all sources for which y is greater than Yo, 
rhether the specific gain is greater than that of the cor- 
esponding uniform distribution or not. Super-gain is 
sed in this wider sense hereafter. 


Yo = (58) 


Si(4mra/d) — 
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The answer to the question of how large a value of 
y/Yo can be tolerated in a particular situation will de- 
pend upon many factors, most of which cannot be con- 
sidered here. It seems reasonable, however, to look upon 
source distributions for which y exceeds Yo by a factor of 
ten or more with extreme caution. The quantity 
(AGo/4a),, is especially sensitive to any shifting of the 
zeros into or toward the visible range and it can undergo 
a ten-fold change with very little provocation. 

Several articles’ !~®° dealing directly or indirectly 
with the subject of super-gain have appeared. Wood- 
ward and Lawson" base their article upon a two-dimen- 
sional formulation of the aperture problem but it con- 
tains concepts closely related to those developed here. 

It remains to discuss F(z, A, %) line sources from the 
point of view of super-gain. If % is less than 2a/, the 
change from the space factor of a uniform distribution 
to an F(z, A, #) space factor involves no migration of 
zeros into or toward the visible range, but only a rear- 
rangement of the zeros already there, and the increase 
in y/Yo will be inconsequential if it occurs at all. Even 
if # is somewhat larger than 2a/\, the relative area 
under the invisible lobes may still be small, especially if 
the design side-lobe level is low. Of course, if % is in- 
creased indefinitely, a point will be reached where 
further increase of # will have little effect upon the 
visible pattern but will continue to raise the level of the 
invisible lobes, making them approach the design level; 
the source can thus acquire an arbitrarily high degree of 
super-gaining. There is, of course, no need to resort to 
large values of # in order to realize most of the beam- 
sharpening potentialities of F(z, A, #) sources, and 
hence there will rarely be any occasion to use such 
sources under conditions which would make y/7» large, 
even though this is theoretically possible. 

By letting % tend to infinity, one approaches the ideal 
space factor, Fo(z, A); the latter, therefore, has an in- 
finite super-gain ratio associated with it. For this reason, 
it is inconceivable that a non-super-gain source could 
out-perform the ideal space factor with regard to the 


11W. W. Hansen, “Notes on Microwaves,” published in limited 
edition by engineers of Sperry Gyroscope Co., Inc., vol. 3, pp. 88- 
100; 1941-1944. 

12S. A. Schelkunoff, “A mathematical theory of linear arrays,” 
Bell Sys. Tech. Jour., vol. 22, pp. 80-107; January, 1943. 

18 C.J. Bouwkamp and N.G. de Bruijn, “The problem of optimum 
antenna current distribution,” Phillips Res. Reports, vol. 1, pp. 135- 
158; 1945-1946. 

14 A. Laemmel, “Source distribution functions for small high gain 
antennas,” Microwave Res. Inst., Polytech. Inst. of Brooklyn, 
Report R-137-47, PIB-88; April 4, 1947. 

16 H. J. Riblet, “Note on the maximum directivity of an antenna,” 
Proc. I.R.E., vol. 36, pp. 620-623; May, 1948. 

16 R, M. Wilmotte, “Note on practical limitations in the directiv- 
ity of antennas,” Proc. I.R.E., vol. 36, pp. 878; July, 1948. 

17 P, M. Woodward and J. D. Lawson, “The theoretical precision 
with which an arbitrary radiation pattern may be obtained from a 
source of finite size,” Jour. IEE, vol. 95, pp. 363-370; 1948. 

18 T, T. Taylor, “A discussion of the maximum directivity of an 
antenna,” Proc. I.R.E., vol. 36, p. 1135; September, 1948. 

19 L. J. Chu, “Physical limitations of omni-directional antennas,” 
J. Appl. Phys., vol. 19, pp. 1163-1175; December, 1948. 

20 N. Yaru, “A note on super-gain antenna arrays,” Proc. I.R.E., 
vol. 39, pp. 1081-1085; September, 1951. 
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compromise between beamwidth and side-lobe level, 
especially if the visible range is extensive enough to in- 
clude at least one complete pair of side lobes. The 
beamwidth-side-lobe relationship given in (51) and il- 
lustrated in Fig. 7 may therefore be regarded as the 
practical optimum compromise for the space factor of a 
line source. The words, “... for the space factor of 

.” are used advisedly, because it is possible to exploit 
the element factor (when there are collinear current 
elements, for example) and use a space factor with an 
increasing side-lobe characteristic to give a flat side- 
lobe characteristic in the total pattern. A small reduc- 
tion of beamwidth can be achieved thus, as Sinclair and 
Cairns?! have shown for arrays of discrete elements, but 
such techniques are beyond this article’s scope. 


THE INVERSE TRANSFORM AND THE 
LIMITING SPECIFIC GAIN 


In calculating the distribution function g(p) to give 
the space factor F(z, A, %), the synthesis method of 
Woodward” will be used. Let g(p) be analyzed into its 
Fourier components: 


fo) 
>> Dnewim? pr Sr? 


g(p) = (60) 


0 Pen 
By applying the Fourier transform (5) to this expres- 


sion, term by term, the following result is obtained: 


Raye toe > pe 1(Z — m) 


m>=—o 


«(z — m) ct 


The terms sin r(z—m)/a(zg—m) are displaced from one 
another by integral steps along the axis of reals in the 
z-plane, hence the principal maximum of any coincides 
with a zero of all others. It follows that F(m) =27D,,, or 


“s F(m) 


iB es 
20 


(62) 


Thus the Fourier coefficients of g(p) are simply the pat- 
tern values at the integers. Since the space factors 
which have been considered are even functions, D_,, 
=D,,, and g(p) in the range p?S7? is given by 


(p) = —{F(0) +25 F(m) cos mpl die onl 62) 
m=1 
If the space factor belongs to the F(z, A, 
is very easy to use because F(m, A, 
m=n, and the summation terminates. 
It is of the utmost importance that the F(m) be cal- 
culated correctly, hence the following simple check may 


n) family, (63) 
n) is zero for 


1G, Sinclair and F, V. Cairns, “Optimum patterns for arrays of 
non-isotropic sources,” Trans. I.R.E., vol. PGAP-1, pp. 50-59; 
February, 1952. 

2 P. M. Woodward, “A method of calculating the field over a 
plane aperture required to produce a given polar diagram,” Jour 
I.E.E., vol. 93, pp. 1554-1558; 1946. 
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be useful. This check is based upon the fact that there 
are two independent methods for calculating g(7), When 
a=0, h(a) =g(a) and, from Theorem III, there is a uni- 
versal asymptotic expression for F(z); 


sin 1Z 
F(z) ~ 2rg(r) as | 2 | > 0, (64) 
TZ 
In other words 

; F(z) 
2rg(r) = lim - (65) 

|z|7 2° SsIn 7z 

TZ 


Substituting the value of 27 g(a) obtained from (63), we 
find that the F(m) must satisfy the following condition: 


F(0) + om 53 (—1)"F(m) = lim te (66) 
sin 72z 


m1 |z| +00 


TZ 


The limiting value of the specific gain is easily ex- 
hibited in terms of the D,,. Appeal is made to the Par- 
seval formula for Fourier series,!° and the iftegral in the 
denominator of (56) becomes 


2 | Dn |?. 


m=—o 


[| ste) bap = 2 (67) 


The integral in the numerator is simply 27D, and the 
gain expression becomes 


(AGo/4a) 53 = 


(68) 


Using this expression, the quantity (AGo/4a)., has been 


plotted as a function of # for four F(z, A, #) line sources 


with design side-lobe ratios of 10, 15, 20, 25 db (Fig. 9). 
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According to (57), when y +1 the limiting value of the 
specific gain is an excellent approximation for the 
specific gain itself. The quantity 7 will be approximately 
equal to unity if the following conditions are simul- 
taneously satisfied: 


1. The ratio 2a/h is of the order of or greater than 
five, for 1<7yoS1.02 when 2a/AZ5. 
2. The source in question has little or no super-gain- 
ing and therefore y=yYo. This will be true for 
F(z, A, 7%) sources when # is of the order of or less 
than 2a/X. 
These conditions roughly define the domain of utility 
of the curves of Fig. 9 from the point of view of estimat- 
ing the specific gain. To be more quantitative, one would 
have to find the indefinite integral of the function 
| F(u, A, n)|*. This hardly seems worthwhile, since the 
majority of practical cases, especially in the field of high- 
gain antennas, will fall well within the domain outlined 
by the two conditions stated above. 


EP 
Lobel 


20 DB DESIGN 
SIDE LOBE RATIO 


As practical examples of the methods proposed in this 
article, complete calculations for 10, 15, 20, and 25 db 
line sources with 7=5 have been performed. The results 
are illustrated in Fig. 10 (preceding page) which shows 
the source distribution and the space factor for each of 
these cases. Fig. 11 shows the corresponding data for a 
20 db side-lobe source with *=10. 
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On the Input Conductance of Thin Antennas* 
GIORGIO BARZILAIt+ 


Summary—Input conductances of thin antennas are computed 
by evaluating separately input voltage and radiated power for a pre- 
scribed value of the current near the maximum. The current distri- 
bution required to evaluate the radiated power and the input voltage 
is obtained by solving Hallén’s equation for the prescribed value of 
the current. First and second order computations were carried out 


for full wave antennas, and the results compared with those ob- 
tained by others. 


INTRODUCTION 


N RECENT years considerable attention has been 
| given to the problem of computing the input ad- 

mittances of thin antennas. With reference to the 
conductance component, two methods of computation 
can be used. The first evaluates the conductance as the 
real part of the ratio of the input current to the input 
voltage. The second computes this quantity as the ratio 
of radiated power to square of applied voltage. 

It is the purpose of this paper to compute the con- 
ductances of thin antennas following the second ap- 
proach, which has not generally been used, and to com- 
pare the results with the existing calculations. 

The evaluation of the radiated power requires knowl- 
edge of the current distribution. This distribution can 
be obtained by prescribing the applied voltage and fol- 
lowing Hallén’s procedure. In our case, however, the 
results are less dependent upon the choice of the ex- 
pansion parameter used to solve the integral equation, if 
the current near the maximum is prescribed. The ap- 
plied voltage must, therefore, be computed. 

First and second order computations have been car- 
ried out for full wave antennas of three different thick- 
nesses. 


THE METHOD OF COMPUTATION 


With reference to Fig. 1, let us consider a thin, center- 
driven, perfectly conducting cylindrical antenna. The 
driving voltage, Vo, is assumed to be applied through 
the gap, d, which is of infinitesimal length. It is well 
known that on the surface of such an antenna scalar 
and vector potentials satisfy transmission line equa- 
tions. These two quantities (if h=/4, where d is the 
wavelength) may be expressed as follows: 


mn 
A(z) = A(h) cos 6 + A(i-—) sin 6 (1) 


oN 
V(z) = jn Ec sin@ — A (i — * cos alee (2) 


where A(z) and V(z) are vector and scalar potentials on 
the surface of the antenna at the height, z. 


* Original manuscript received by the PGAP, March 31, 1954; 
revised manuscript received October 11, 1954. 
+ Polytechnic Inst. of Brooklyn, Brooklyn, INE Ye 


n = V(u/e) & 1207 in free space, 


Mande are the permeability and the dielectric constants 
of the medium surrounding the antenna, 
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Fig. 1 


On the other hand, the vector potential can be ex- 
pressed in terms of its sources, the conduction current 
densities. In the case under discussion we can approxi- 
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mate A(z) as follows 


jkr 


1 +h e 
A(z) = — { 
\ Ar —h if 


r= [a+ (2 — 2’)? ]/, 


where J(z) is the current in the antenna at the height, z. 

Upon taking (3) into account, (1) becomes an integral 
equation for the current which can be formally solved 
following Hallén’s procedure by letting 


I(z’)dz’ (3) 


1 
A(e) = = [2@)2 + fle) (4) 
2h 
Q = 2in — (5) 
a 


where f(z) is a function essentially independent of a. 
We obtain 


\ 
I(z) = 1(1 = *) sin 6 
= =| cos 6 — f(z) AG = *) sin | . (6) 
Q ; 4 


Assuming now I|h—(A/4)] prescribed, Z(z) can be ex- 
panded in a series of descending powers of 2 by using 
the well known iteration procedure.! As zero order solu- 
tion we shall take 


Ih = 2 i 
o(Z) = 1(1 — *) sin 6. Gp) 


It has to be noted that solutions of any order vanish 
identically for z= +h. 

The first order solution, immediately obtained from 
(6), is 


T(z) =Io(z) 
mene | h) cos é— + pes i | 8 
fi(h) cos 8—fi(z) i ) sin Ree oy) 


where f(z) is obtained from (4) assuming the current 
distribution, (7). Higher order solutions can be written 
down readily if desired. 

We now recall that 


Vo= 2V(0"), (9) 


where V(O*) is the value approached by V(z) when 
z—O from the positive side. From (2) we obtain 


oN 
Vo = 29 E (h) sinkh — A (i — *) cos an. (10) 


By using (7) for computing A(h) and A[h—(\/4)] a 
first order solution, Vo, is obtained for the applied 
voltage, Vo. By using (8) a second order solution, Vo, 
can be obtained. 


‘ E. Hallén, “Theoretical investigations into the transmitting and 
receiving qualities of antennas,” Nova Acta (Upsala, Sweden), vol. 
11, no. 4; 1938. 
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We notice that since V(z) can be expressed in terms 
of constants other than A(h) and A[h—(h/4)], other 
expressions for V» can be obtained. For instance we can 


write 
0 if) . 


If a certain order current is used to compute Vo, (10) 
and (11) will not, in general, yield same value for Vo. 

The power radiated can be computed by the following 
expression 


(11) 


+kh +kh 
K(z — 2’)I(z)I*(2’)dzdz’ 


—kh —kh 


P = 30 (12) 


where the asterisk indicates the complex conjugate and 
1 /sin u 
VG OA ee GOS |e 
ur\ u 


Eq. (12) lends itself to computations via the method of 
moments to yield approximate results. A detailed ac- 
count of this method, as well as a justification of (12), is 
given by Schelkunoff and Friis.? 

Once V, and P are determined the input conductance, 
G, is immediately obtained by 


2 
VoVor 


G= (13) 
NUMERICAL COMPUTATIONS 


All the numerical computations have been carried out 
for h=)/2. I(h—X4) =1 has been assumed. Formulas to 
compute 


Ala) = fil(z) + gf" (2) 
are given in the literature.* We found 


fila) = — 3.5523 — 71.9405. (14) 


The first order voltage, Vo, is obtained from (10) or 
(11) by taking into account (4): 


Vor = 60[1.9405 — 7(Q — 3.5523) ]. (15) 
The numerical values for the first order powers are given 
in the literature.‘ In our case we have 


P, = 99.54. (16) 


The first order conductance, Gi, can therefore be ex- 
pressed as follows 


2 199.09 
G, = z = i (17) 
K,.(Ka — 186.28) + 22231 


2S. A. Schelkunoff and H. T. Friis, “Antennas—Theory and 
Practice,” John Wiley and Sons, Inc., New York, N. Y., 1952. 

* R. King and D. Middleton, “The cylindrical antenna: current 
ane ved ance Quart. Appl. Math., vol. 3, pp. 302-335; January, 

*S. A. Schelkunoff, “Advanced Antenna Theory,” John Wiley 
and Sons, Inc., New York, N. Y., 1952. 
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In (17) G, has been expressed in terms of the average 
characteristic impedance for cylindrical antennas de- 
fined as follows: 


K, = 60(0 — 2). 


to facilitate the comparison with the corresponding 
formula obtained by Schelkunoff, using his mode theory 
of the antenna.‘ 

To carry out second order computations, the second 
order voltage, Vos, and the second order power, Po, 
have to be computed. The first order current necessary 
to compute these quantities is obtained from (8). It is 
interesting to notice, at this point, that the difference, 
I,(z) —Io(z), which enters into all the calculations of the 
second order quantities, is inversely proportional to Q. 
Furthermore we notice that if, instead of Q, we use an 
expansion parameter, Y = Q+H, where H is an arbitrary 
constant, the quantity inside the brackets in (8) will not 
change its value. This is apparent if we substitute Q=y 
— fH into (4), and we recall that (8) has been obtained 
using the current distribution, (7). These characteristic 
features of the method under discussion simplify the 
computations of Vo. and P» for different values of Q and 
different expansion parameters. 

In what follows Q will be assumed as the expansion 
parameter; however, all the results obtained will still be 
valid if w is substituted for Q. 

Voz can be computed from (3) and (11), and for the 
desired approximation simple numerical integrations 
are quite adequate. 

The sinusoidal and cosinusoidal components appear- 
ing in (8) can be integrated without difficulty. To com- 
pute the contribution of the fi(z) component to V2, we 
have applied the trapezoidal rule to the function, 


PO 
Ale) - 4 (>) 


Two sets of computations have been carried out, corre- 
sponding to the division of the interval of integration 
into 16 and 32 equal parts. 

For 16 parts we obtained 


C= jkr 


- 
- 


60 
Voo — Va = ee (Se ee» (18) 
For 32 parts we obtained 
60 


The error made by using (19) has been estimated to a 
few units in the last decimal place of the quantities in- 
side parenthesis. This approximation is sufficient for our 
purposes as will be discussed later. 

The computations for the second order power have 
been carried out by using the method of moments. This 
method consists essentially of replacing the integral, 
(12), with a double sum. It is equivalent to dividing the 
‘antenna into a certain number of sections and replacing 
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each section with an elementary dipole having a mo- 
ment equal to the integral of the current along the sec- 
tion under consideration and situated at the center of 
the same section. 

It isa simple matter to integrate f;(z) in (8) analyti- 
cally. Once the integral of f(z) is known, the moments 
required for further computations follow immediately. 
For h=X/2 we obtained 


2k i fi(z) dz 
0 


= [—2 Si 2kz+Si 2(r— kz) +Si 2(r+ kz) | sin kz 
+ |—2 Ci 2kz— Ci 2(r— kz) + Ci 2(r+ kz) | cos kz 


TZ N/2+2 
= |- ay Jina ln | cos kz 
IN d/2 


ae 


+4 Ci kz—2 Ci (m—kz)+2 Ci (r+hz) 
+7} [—2 Ci 2ke+Ci 2(r—kz) + Ci 2(r+hz) | sin kz 
+ [+2 Si 2kz+Si 2(r—kz) —Si 2(r+hz) | cos kz 


oes 
+ | +n (—-1)| sin kg + 


—4 Si ke +2 Si 2(r—kz)—2 Si 2(n+hz) } 


where y=5//2... and 220. 

Two sets of computations have been carried out for 
the second order power corresponding to the division 
of the antenna into 4 and 8 equal parts. 

With 4 divisions we obtained 


(20) 


1 1 
P= Py ae eer eet (21) 
OF Q 
P, = 98.05. (22) 
With 8 divisions we obtained 
1 1 
1 aang es OS Ci) (23) 
2 Q 
P, = 99.12. (24) 


By comparing (16), (22) and (24), and noting that (21) 
and (23) yield values which differ by about 5 per cent 
for practically any value of Q, we have estimated the 
error made in computing P:—P, by (23) to be a few per 
cent. We shall see that in this case also this approxima- 
tion is sufficient for our purposes. 


TABLE I 
Q Ge Gio GoW Got Go GS 
10 1220 1.081 12033 0.869 0.995 1.243 
15 0.410 0.383 0.376 0.345 0.374 0.403 
20 0.202 0.192 0.190 0.181 0.190 0.197 


In Table I we have recorded the results of our com- 
putations together with a list of the symbols used. In 
all the computations we have assumed for P; the value 
given by (16). King and Middleton's expansion param- 
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eter, W for h=X/4, is defined as follows: 


= («-=)| 
y Saget h adhe 


As far as the accuracy for Ge,g and Gy is concerned, 
we notice that if (18) instead of (19) is used, the last 
decimal place of G2, will change by one unit for Q=10, 
and will undergo no change for the other two values of 
Q. We notice also that if (21) instead of (23) is used, the 
last decimal place of G2, will change by three units for 
Q=10, by one unit for Q=15, and will undergo no 
change for 2=20. It therefore seems reasonable to con- 
clude that for Q=10 an error of two or three units may 
affect the last decimal figure of G2,9 and G2,y while for 
the other two values of Q the figures recorded in Table 
I should be correct. 

Inspection of Table I shows that the method discussed 
leads, at least as far as second order results are con- 
cerned, to values for the input conductance which are 
rather insensitive to the choice of the expansion parame- 
ter. Go, and», Gy do, in fact, differ by about four per cent. 
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Conductances are in millimhos. The symbols, G, rep- 
resent conductances computed as ratios of radiated 
power to the square of the input voltage. The symbols, 
G, represent conductances computed as the real part of 
the ratio of input current to input voltage. 


List OF SYMBOLS 


G,=First order conductance computed by 
using (17). 

Cee and G:=Second order conductances computed by 
the method discussed in this paper, using 
the expansion parameter Q and w respec- 
tively. 

G,¥ and G,X =Second order conductances according to 
Hallén,® King,*® and Middleton.® 
GS =Conductance according to Schelkunoff’s 
mode theory. 

'E, Hallén, “Admittance diagrams for antennas and the relation 
between antenna theories,” Cruft Laboratory, Harvard University 
Report 46, June 1, 1948. : 

6p, A. Kennedy and R. King, “Experimental and theoretical 


impedances and admittances of center-driven antennas,” Cruft Lab- 
oratory, Harvard University Report 155, April 1, 1953. 


Theory of Radio Reflections from Electron-Ion Clouds* 
Von R. ESHLEMANT 


Summary—Approximations to the reflection coefficients of elec- 
tron-ion clouds of various sizes, shapes, and densities are determined. 
A simplified approach to the problem is employed, wherein the total 
reflection from the low-density clouds is computed by summing the 
scattered wavelets from the individual electrons, while the high- 
density clouds are treated as total reflectors at the critical density 
radius. Some of the limitations of this method are discussed. While 
more elaborate determinations should be used in certain regions of 
cloud size and density, the method used here provides an over-all, 
first-order approximation of the effects of size, shape, and density 
on the reflecting properties of electron-ion clouds. Several possible 
applications of the theoretical results are outlined. 


INTRODUCTION 
A THEORETICAL investigation of radio reflec- 


tions from ionized regions has direct application 

to the study of such ionospheric irregularities as 
sporadic-# ionization, auroral ionization, and meteor 
ionization trails. An attempt is made here to obtain 
results which are general enough to provide a first-order 
approximation of the effects of size, shape, and density 
on the reflecting properties of electron clouds. More 
thorough analyses have been presented elsewhere for 


* Original manuscript received by the PGAP, March 31, 1954; 
revised manuscript received, October 11, 1954. This work was sup- 
ported jointly by the U. S. Navy (Office of Naval Research), the 
U.S. Army Signal Corps, and the U. S. Air Force, Contract N6onr- 
ole asians 

} Radio Propagation Lab., Stanford Univ., Stanford, Calif. 


ionized regions having more specific properties. How- 
ever, it is believed that the simplified treatment which 
follows may result in a useful method for analyzing the 
salient characteristics of radio reflections from electron- 
ion clouds having general sizes, shapes, and densities. 
Only the back-reflected signals are determined; that is, 
the radio transmitter and receiver are assumed to be at 
the same location. 

In order to determine the general effects of cloud 
shape, the three basic geometrical configurations are 
considered. That is, the electrons are distributed in 
dimensions about a (3—m) dimensional origin. For 
n=1, n=2, or n=3, the electron volume density is a 
function of the distance from a reference plane, line, or 
point, and the clouds have a laminal, cylindrical, or 
spherical shape, respectively. A simplified method is 
employed to determine the reflection from such clouds 
having any size and electron content, and having several 
different distributions of electron density as a function 
of distance from the origin. Particular attention is 
given the clouds which result from normal diffusion in 
n dimensions from an instantaneous source of (3—2) di- 
mensions. 

It is assumed that for low electron densities, each 
electron scatters freely and independently, while for 
high densities, total reflection occurs at the surface of 
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critical electron density. The reliability of these as- 
sumptions is discussed. The effects of the geomagnetic 
field, recombination and attachment, and collision and 
radiation damping are not considered. These factors are 
usually unimportant for high-frequency radio reflec- 
tions from ionospheric irregularities in the E-region. 

The emphasis in this paper is on the development of 
the theoretical results, though several possible applica- 
tions are suggested. 


Low-DeEnsity CLoups 


The scattered electric field strength at range R’ from 
a single electron moving in the £ direction with accelera- 
tion £ is, in meter-kilogram-second units, poet sin 
a/4mR’, where io is the permeability of free space, é is 
the magnitude of the electronic charge, and a is the 
angle between the é direction and the direction to the 
field point. The acceleration of a free electron acted 
upon by an incident electromagnetic wave is, omitting 
the exp (—iwt) time variation, —Ein.e/m, where m is 
the electron’s rest mass. In the problems under con- 
sideration, the radio transmitter and receiver are at the 
same location, so that sin a=1. Thus, the ratio of the 
scattered-field strength at the receiver to the incident- 
field strength at the single electron is given by 


Be Loe? 1 
= -(~)—., (1) 
4m / 7R’ 


eine 
where poe?/4m =0.8852-10-" meters. With the assump- 
tion that the electrons scatter freely and independently, 
the scattered-field strength from an ensemble of elec- 
trons may be found by vectorially adding their indi- 
vidual contributions. When the extent of the cloud is 
small compared to its range, the incident wave is essen- 
tially constant in amplitude throughout the cloud. Thus, 


Pre (; =) : fp N exp (2ikR’)d (2) 
=> —_ SS SS That ©. 4 4 vy, 
Pre 4m wR » e 


where R is the range to a reference point within the 
cloud, N is the electron volume density, and & is the 
wave number, 27/2. 

In Fig. 1, laminal, cylindrical, and spherical electron 
clouds are illustrated. While these clouds are infinite in 
some dimensions, the phase factor controls the echo 
amplitude in such a manner that the important scatter- 
ing region is small in extent, so that (2) still applies. 
Note that the co-ordinates are chosen and named so 
that, in each case, N is a function of 7 only. Assuming 
that r<(mR/k)*”, (2) can be simplified to read 


ola) | 
Eine z 4m - 2 
N(r)F,(r)r*—1dr 


-28-2(gkR) I”) (2 - (3) 
f N(v)r2 dr 
0 
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In this expression, g, is the total number of electrons 
/meter>” of the cloud and F,(r) =(kr)@-”?D(n/2) 
J (n-2)j2(2kr), so that F,(r) =cos (2kr), F,(r) = Jo(2kr), 
and F;3(r) = [sin (2kr) |/(2kr). 

Let the dimensionless quantity which is bracketed in 
(3) be represented by X,. Thus, Xi is proportional to 
the wavelength and to the number of electrons per 
square meter in the laminal clouds. Similarly, X»2 is 
proportional to the number of electrons per meter in 
the cylindrical clouds, and X3 is proportional to the 
number of electrons in the spherical clouds and in- 
versely proportional to the wavelength. 


n=] 
laminal 


n=3 
spherical 


Fig. 1—Laminal, cylindrical, and spherical electron clouds. 


It is convenient to define the reflection coefficients pn 
of the electron clouds in the following manner: 


Esc 
—* = 2-8) 12 kR/2) O12, (4) 


ine 

The (rkR/2)°-”” factor is a normalizing expression 
commonly used to make the geometrical optics reflec- 
tion coefficients of metal cylinders and spheres inde- 
pendent of range and proportional to some function of 
the target size in wavelengths. The 2‘"-®/? factor orig- 
inates from the difference between plane-wave and 
spherical-wave incidence upon the target. Other ex- 
pressions which are often used to describe the reflection 
properties of radar targets may be related to the reflec- 
tion coefficient defined above, from the following defini- 
tions: (E,./E ine)? =(2RR/G)? Pr/Pr=o/47rR?, where G 
is the antenna power gain over an isotropic radiator, 
Pr and Pr are the received and transmitted peak 
powers, and g@ is the radar cross section. 

Using (4), the reflection coefficients of the low-density 
clouds become 
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TABERMI 
REFLECTION COEFFICIENTS FOR Low-DENsITY CLOUDS 


Radial electron 


. pa 
distribution Puls pr 3 
Homogeneous 
N/No=1, 0<r<a (X1 sin b)/b 2XoJi(b) /b 3X;3(sin b—b cos b)/b8 h 
=0,a<r< [o-] [b-3/2] [0] 
Linear 


N/No=1—r1r/a, 0<r<a 2X,(1—cos b)/b? 


[o-?] 


=0,a<r<o 


6x} 100) 42/0) Se Som s1(b) ¢ /0? 


12X3(2—2 cos b—b sin b) /b* 


[o-*] 


[o-o 


Parabolic 
N/No=1-—(r/a)?, 0<r<a 


3Xi(sin b—b cos b)/b? 
=0,a<r<e% [b-?] 


8X2Jo(b) /b? 
[a5] 


15.X3(3 sin b—30 cos ch sin 5) /b® 
6-3 


Exponential 


N/No=exp (—r/a) Cy eae 


X3/(1+6?)? 


Xo/(1+02)3/ 
[o-8] [>-*] 


Gaussian 
N/No=exp (—71?/a?) X, exp (—5?/4) 


[exp (—0?/4) | 


X_2 exp (—b?/4 


) X; exp (—b?/4) 
[exp (—b?/4)] 


[exp (—b?/4)] 


In these expressions, b=2ka and Xn = (joe?/4m)qn(k/2)"-*; where k=27/, (ue?/4m) =0.8852-10~'4 m, a=cloud size as defined in the ; 
expressions for the radial electron distribution, gi =electrons/m? in the laminal cloud, g:=electrons/m in the cylindrical cloud, gs =electrons 
in the spherical cloud. The quantities in brackets indicate the dependence of paz on 6 for large b. 


f NOE AD dr 


Hl N(r)r"— dr 
0 


If the extent of the cloud is small in the 7 direction, 
Pnt=X, regardless of the functional behavior of NV(r) 
since F,,(r)&1 for rn. 

Eq. (5) has been solved for five different NV functions, 
and the results are tabulated in Table I. The results 
for n=2 have previously been presented by several 
authors.!- In the table, = 2ka, and a is the cloud radius 
as defined in the expressions for NV. Note that the max- 
ima of pnz are equal to X, and occur when )=0. For the 
homogeneous, linear, and parabolic models, the func- 
tional relationships between the reflection coefficients 
and the sizes of the clouds exhibit nulls and lobes. The 
envelopes of these patterns decrease monotonically 
with increasing cloud size in the manner given by the 
functions in brackets which accompany the expressions 
for pnz in Table I. The reflection coefficients of the ex- 
ponential and Gaussian models decrease smoothly with 
increasing cloud size, and their final dependence on size 
is given by the expressions in brackets. 

The function V may be expected to vary with time 
due to possible time variation of the electron sources 
and sinks, and due to diffusion, wind shear, and tur- 
bulence. While the reflection coefficient may be found 
from (5) for any given JN, the time behavior of NV and 
hence pnz is, in general, very difficult to determine. For 


(5) 


PnL = Men 


1]. Feinstein, “The interpretation of radar echoes from meteor 
trails,” Jour. Geophys. Res., vol. 56, pp. 37-51; March, 1951. 

2, N. Herlofson, “Plasma resonance in ionospheric irregularities, ” 
Arkiv for Fysik, vol. 3, pp. 247-299; April, 1951. 

3 V. R. Eshleman, “The Mechanism of Radio Reflections from 
Meteoric Ionization,” Tech. Rep. No. 49, Electronics Res. Lab., 
Stanford Univ., Stanford, California; July, 1952. 


the simple case of an instantaneous (3—m) dimensional 
source and no electron sinks (i.e., no recombination or 


: 
| 


attachment), the forces of diffusion cause the radial dis- 


tribution of electrons to be Gaussian as long as there is 
no disturbing wind shear or turbulence. Notwithstand- 
ing the restrictions in this picture of cloud formation 
and growth, many of the characteristics of the echoes 
received from meteor ionization trails (where =2) 
have been explained by use of the Gaussian model.*-* 

For the reasons cited above, the clouds which result 
from uniform expansion from an instantaneous plane, 
line, or point source are of particular interest. For such 
clouds, the normal diffusion equation is 


oN 0 oN 
= De = (= —), (6) © 
at or or 
with solutions’ 
V(r, 0) a ( Ss ) (7) 
N(r, t) = ———— exp | — —), 
(4nDi)2 * \~ api 


where D is the diffusion coefficient in square meters per 
second and f is the time in seconds from the instant of 
formation of the source. Thus, the electron volume 
density has a Gaussian distribution about the source, 
and, from Table I, 


(pnt)g = Xp exp (— 4k2D2). 


‘T. R. Kaiser and R. L. Closs, “Theory of radio reflections from 
meteor trails: I,” Phil. Mag., vol. 43, pp. 1-32; January, 1952. 
5 J. S. Greenhow, “Characteristics of radio echoes from meteor 


trails: III. The behavior of the electron trails after formation,” Proc. 


Phys. Soc., vol. 65, pp. 169-181; March, 1952. 


(8) | 


q 


et Kaiser, “Radio echo studies of meteor ionization,” Re-_ 
search Reviews (Phil. Mag. Supplement), vol. 2, pp. 495-544; October, — 


1953. 


7H. S. Carslaw, “Introduction to the Mathematical Theory of — 


Conduction of Heat in Solids,” Dover Press, 2nd ed., pp. 149-153; 
New York, N. Y. 
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TABLE II 
REFLECTION COEFFICIENTS FOR H1GH-DENSITY CLoups 


Radial Electron 
Distribution PlH 


p2H Pat 
aN. Homogeneous | 
o=1,0<r<a 1 (xb /8) 1/2 
ag ¥ f b/8) 
0,a<r<e X1>b/2 X2>1b?/16 eu X3>71b3/48 
Linear 
V/No=1—r/a, 0<r<a 1 (xb 2 2 71] 
be FY /8)"/?(1 —b?/48X2)1/? b/8) (1 —3b3/192X 
=0,a<r<o X,>6/4 Xo >1b?/48 ey Dy gee lores | 
Parabolic 
ez) O0<r<a ete _ (wb/8)1/2(1 —1b2/32X2) 1/4 (xb /8)(1 —2b3/120X3) V2 
SES ee X1>b/3 X_>nb?/32 X3>7b/120 
Exponential 
V/No=exp (—r/a) i (3b/8)!/2[In (8X o/b?) ]/2 (1b/8) In (8X3/1b*) 
Xi>6/2 X_>71b?/8 X3>7b3/8 
M/N Gaussian 
o=exp (—?r?/a?) il (xb/8)2/2[In (16X>/ab2) ]1/4 b/8) 11 3/23 |1/2 
Xy > 71!2h/4 pe | TN Sete Ae! 


X3> 13/23/64 


_ In these expressions, b=2ka and Xy, = (uoe?/4m)qn(R/2)"-2; where k=27/d, (uoe?/4m) =0.8852 - 10-147, a=cloud size as defined in the ex- 


sressions for the radial electron distrib 
che spherical cloud. 


ution, g:=electrons/m? in the laminal cloud, g:=electrons/m in the cylindrical cloud, g;=electrons in 


The expressions at the bottom of each box indicate the values of X, needed for a critical-density radius to exist in the clouds. 


The subscript G on pyz indicates that the N function is 
Gaussian. Note that the maximum values of the reflec- 
tion coefficients occur at zero time. The echoes then de- 
cay exponentially with a time constant, given by (Tn1) ¢ 
= (4k?D)—, which is independent of the electron density 
and the configuration of the cloud. 


H1GH-DENsITy CLoupDSs 


It is well known that an electromagnetic wave cannot 
propagate in an ionized gas when the electron density is 
so large that the equivalent dielectric constant is nega- 
tive. If the region of negative dielectric constant is large 
enough in extent for the evanescent (nonpropagating) 
Wave to be appreciably reduced in amplitude, an ap- 
proximation to the radio reflection from the electron 
cloud is obtained by assuming that total reflection 
occurs at the surface where the equivalent dielectric 
constant is zero. 

The critical volume density (the electron density cor- 
responding to an equivalent dielectric constant of zero) 
is given by N.=(k/2)?/(uce?/4m). The simultaneous 
solution of this equation with the NV equation for the 
electron cloud gives the critical density radius of the 
high-density cloud. Assuming that the geometrical op- 
tics equations for reflections from plane, cylindrical, and 
spherical mirrors apply to this problem, the scattered 
field from a perfect reflector at the critical density ra- 
dius fre is given by Eec/Eine=(fne/R)” ??/2. The 
validity of this assumption will be discussed in a later 
section. 

The reflection coefficients pnrx of the high-density 
clouds having the same WV functions as were tabulated 
in Table I for the low-density clouds are presented in 
Table II. With each value of pry is given the value of 
X,, required for a critical density radius to exist in the 


cloud. That is, a criterion is presented for determining 
whether a cloud is of high or low density. It will be 
shown later that there is an additional requirement for 
this determination which must be used in conjunction 
with the expressions given in Table II. 

For the expanding Gaussian clouds with high electron 
density, the time dependence of the reflection co- 
efficients may be found from (7) and Table II. Thus, 


eee fee oe: amie ‘ 
Pr) Gg = A n lear : ( ) 


As the cloud expands, the critical density radius in- 
creases to a maximum value and then decreases to zero. 
The time to the disappearance of the critical density 
radius is given by (Ina) g=X7,7!"/rk?D. After this time, 
the low-density formula should be used to find the re- 
flection coefficients. However, it will be shown that in 
addition to the values given in Table II, X, should be 
greater than about unity in order for the high-density 
formulas to apply. With X,>1 and t>(Tirz)e@, (par)e 
is vanishingly small, so that (T,#) q@ can be taken as the 
total duration of the echoes received from the high- 
density Gaussian clouds. 


MAXIMUM REFLECTION COEFFICIENTS 


The maximum values of the reflection coefficients with 
constant X,, and the cloud sizes at which these maxi- 
mum values occur, are of some interest. It has already 
been shown that for any N function the maximum value 
of pnt is Xp, and that this value occurs at zero size. In 
Table III on the following page, the maximum values of 
pnu are listed. Also listed in Table III are the values of b 
at which these values of reflection coefficient occur. 
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TABLE III 
MaximuM REFLECTION COEFFICIENTS OF H1GH-DENSITY CLOUDS 
| 
oe cues (pix) max (p2H) max (ps1) max 
Homogeneous 
AN Oy 1 (wX2/4)1!4 (322X3/32)1/3 
(0), Sirk E29 b6< 2X1 b?=16X2/a 63 =48.X3/r 
ree | 
Linear te 
N/No=1-r/a, 0<r<a 1 (1X2/9)1/4 (97?2X3/128) | 
=), <i K 9 b< 4X, 62?=16X9/r b3=48X3/r : 
Parabolic | 
N/No=1—(r/a)?, 0<r<a 1 (4 X/8)1/4 [(3/5)1/2942X3/160] 1/8 
==) <p <0 b<3X, b?=16X0/r 63=48X3/mr 
Exponential 
N/No=exp (—r/a) 1 [rX>/2 exp (2) ]1/4 [2772X3/64 exp (3) ]*/8 
b<2X, b?=8X.2/zm exp (2) b3 =8X3/m exp (3) 
Gaussian is 
N/No=exp (—1?/a?) 1 (rX2/4 exp)!/4 [(3x/2)3/2X3/8 exp (3/2) }¥/ 
b<4X,/ri!? b?=16Xe/7 exp b3 =64.X3/73/2 exp (3/2) 


In these expressions, b=2ka and Xn=(woe?/4m)qn(k/2)"-?; where k=2m/X, (uoe?/4m) =0.8852-10-4 m, a=cloud size as defined in the 
expressions for the radial electron distribution, g:=electrons/m? in the laminal cloud, g2=electrons/m in the cylindrical cloud, qs =electrons 


in the spherical cloud. 


The expressions at the bottom of each box give the cloud size at which the maximum reflections occur. 


In Table IV (below), wavelength, range, and elec- 
tron-density dependencies of maximum reflection co- 
efficients are summarized. The same dependencies of the 
maximum values of E../Einc, Pr/Pr, and o are also in- 
cluded for the convenience of those who are better 
acquainted with one of the alternative expressions for 
the reflection properties of the target. As can be seen 
from Table III, the functional dependence of the max- 
imum reflection on the wavelength, range, and electron 
density is independent of the radial electron distribution 
in the cloud. 


DIvis1on BETWEEN HIGH- AND Low-DENSITY 
CLOUDS 


In Fig. 2, the square of the maximum reflection co- 
efficients obtained from Tables I and III are plotted as 
functions of the dimensionless parameter X,. The single 
line for (pn1r)*max intersects the lines for (pnz)’max Near 
X,=1. This result, together with exact analyses for 
several cloud models, indicates that X,=1 marks the 
lower limit of applicability of the high-density formu- 
las. That is, with X,<1, the low-density formulas 
should be used regardless of cloud size. While a critical 


Fig. 2—Maximum reflection coefficients plotted as a function of Xn, 


where Xn = (uoe?/4m) gn (k/2)"-*. Positions within the shaded 
regions depend upon the cloud model (see Table III). 


TABLE IV 
MaximuM Ecuo DEPENDENCE ON }, R, AND Gn 
Fs Cloud Low density High density 
configuration p Ecc/ ine Pr/Pr o p Esc/ Eine Pr/Pr o 
5 Laminal AIR%qyt ARG MR~q)? MR2q? ORG? PRG? Rg)? PRAM? 
wees Cylindrical NPRGs! AU/2R-1/ 201 MR-3qo? AIR" go? NOR%21/4 MARU /4 | AIR—3get/2 Rigel? 
3 | Spherical ATER%|Qs! | ALR1g51 MWR-4g3? ORGS? ATN3R0ga1/3 d2/3R—-19,1/8 A10/8R-4952/8 A4/3R0g,2/8 


Se st 


Maximum echo dependence on }, R, and qn; \=wavelength, R=range, qi =electrons/m? 


cylindrical cloud, and g;=electrons in the spherical cloud. 


in the laminal clouds, g:=electrons/m in the 
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lensity may exist in such a cloud, the region of negative 
lielectric constant is not large enough for the incident 
vave to be appreciably attenuated. With X,, greater 
han unity but less than the values given in Table MM; 
he low-density formulas should be used, since there are 
10 critical density surfaces in the clouds. The high- 
lensity formulas apply, therefore, only when X,, is 
yreater than unity and also greater than the values 
isted in Table II. 

While a sharp break between the low- and high- 
lensity regions has been shown in the curves of Pigx2; 
t is expected that a more detailed treatment would in- 
licate a gradual transition between the two regions. 
For the homogeneous slab (n=1), the exact treatment 
zives the maximum value of p; as tanh X. This function 
joins the low- and high-density regions in a smooth man- 
ner. For n=2, exact solutions have been obtained for a 
nomogeneous cylinder,*’ and they again show a smooth 
transition between the two lines in Fig. 2 which apply 
to the homogeneous cylinder. 


Fig. 3—Reflection coefficients of homogeneous, laminal clouds by 
wave matching (---) and approximate (—) analyses. The exact 
(wave matching) curves are plotted from 


De b? — 20X, ba 
aw Xi? sin? [(62 — 26X1)"2]) 
For X,=0.5 (low-density), the approximate formula is 


pit =(X1 sin b)/b, while for X:=5.0 (high-density), maa=1, 
0<b/2< Xi; pin =(X1 sin b)/d, Xe 0/20, 


To illustrate the use of the approximate formulas for 
X, less and greater than unity, the reflection co- 
sfficients for a homogeneous slab with X,=0.5 and 
¥,=5.0 are plotted in Fig. 3. Both the exact and ap- 
proximate solutions are given for comparison. 


8V. R. Eshleman, “The effect of radar wavelength on meteor 
cho rate,” Trans. I.R.E., vol. AP-1, pp. 37-42; October, 1953. 
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Ecuo DURATIONS 


For the Gaussian clouds expanding by normal dif- 
fusion, the time behavior of the echoes is given in (8) 
and (9). Following these equations are expressions for 
the echo durations. In Fig. 4, normalized expressions for 


loop 


is 
4kD(Tg 


Xn : 


=i ie 
107! ! 10! 10 


10 


Fig. 4—Echo durations of Gaussian clouds. (Tnr)¢=1/4k2D and 
(Tra) @=X12!/"/rk2D, where D is the diffusion coefficient and 
Xn = (moe?/4m) gn (k/2)". 


the echo durations are plotted as a function of the di- 
mensionless parameter, X,. The horizontal line applies 
for all m in (Tnz)¢. The three sloping lines are plots of 
(Tnu)e¢ for the three possible values of 1. These three 
lines intersect the horizontal line near X,=1. More de- 
tailed computations in this region should show a gradual 
transition between the curves for the two density re- 
gions, as well as a smooth transition between the echo 
shapes represented by (8) and (9). 

All the echo durations for Gaussian clouds are in- 
versely proportional to the diffusion coefficient, and also 
are dependent upon wavelength and q,, as listed in Ta- 
ble V (next page). Only Gaussian clouds maintain their 
identity under the forces of diffusive expansion. How- 
ever, Table V would also apply to the other cloud 
models if there were some other expansion mechanism 
which would change only 0 in the V/N functions, and if 
b? were proportional to ¢ as it is for the Gaussian clouds. 


LIMITATIONS OF THE SIMPLIFIED APPROACH 


It is desirable to obtain some indication of the reli- 
ability of the method which was employed above to 
find the reflection coefficients of various electron clouds. 
In the low-density region, more detailed analyses'* 
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have shown that under certain conditions, polarization 
forces may be set up within the clouds. That is, the 
motion of the electrons in a direction of changing den- 
sity sets up electric fields which act as restoring forces 
on the electrons. For n=1, and for n=2, if the electric 
vector is parallel to the axis of the cylindrical cloud, 
there can be no polarization forces since there is no 
density gradient in the direction of electron motion. 
However, for n =3, and for n =2, with the electric vector 
perpendicular to the axis of the cylindrical cloud, the 
electron motions may be no longer mutually inde- 
pendent. With a restoring force on the electrons, they 
have a natural period of vibration. If the incident radi- 
ation has a period corresponding to a weighted mean of 
the natural periods of the electrons, the reflection co- 
efficient will be greater than that computed by assum- 
ing mutually independent scattering. This effect has 
been called “plasma resonance” by Herlofson.? The 
magnitude of the resonant response is dependent upon 
the diffusiveness of the cloud, the response being 
greater for the clouds with the larger density gradients 
at the critical density radius. The homogeneous clouds 
have the largest possible density gradient at 7n-, and 
the reflection coefficient here may be several powers 
of 10 greater than that computed by the method 
given above. For the Gaussian clouds, however, Kaiser 
and Closs* have shown that the maximum resonant 
response is only 2 for n=2, and is essentially non- 
existent for n=3. In the low-density clouds, the 
polarization forces are no longer important if b>2, and 
all the resonant responses occur at smaller cloud sizes. 
For cloud radii less than the resonant sizes, the re- 
flections are less than those computed from independent 
electron scattering, since the same binding forces which 
cause resonance now restrict the motions of the elec- 
trons to less than they would be under the influence of 
the incident wave only. 


TABULA, 


Ecuo DuRATION DEPENDENCE ON \ AND fis 
FOR THE GAUSSIAN CLOUDS 


Low density High density 
n (duration to some percentage | (total duration of critical 
of initial amplitude) density) 
1 qu \Mq.2 
2 d?q2° qa" 
2 2939 \4/3g,2/3 


X\=wavelength, q:=electrons/m? in the laminal cloud, g2=elec- 
trons/m in the cylindrical cloud, and g;=electrons in the spherical 
cloud. 

In the high-density clouds, the resonant response is 
negligible since, at the size required for resonance, the 
polarization forces are out of phase across the depth of 
the cloud. In addition, there are other possible sources 
of error in the simplifying assumptions which were made 
for the high-density clouds. In effect, all of the ionization 
outside the critical density radius was replaced by a per- 
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fect reflector at this radius. However, Manning® has 
shown that the ionization beyond the critical density 
radius defocuses the radiation. The magnitude of this 
effect depends upon the diffusiveness of the cloud. There 
is no defocussing for the homogeneous model since there 
are no electrons outside the critical density radius. For 
the Gaussian models, the maximum reflection coefh- — 
cients given in Table III should be multiplied by the 
factor 2¢-”/2 to account for this effect. 

In the preceding analysis, a geometrical optics ap- 
proximation has been used. It should be pointed out 
that this approximation is not valid for perfect reflectors 
of small size. With b<2, the reflection from a metal 
sphere is given by the Rayleigh Law which departs 
rapidly from the geometrical optics value as b de- 
creases. Also, with ) <2, the perpendicular and parallel 
reflection coefficients of a metal cylinder differ from the 
geometrical optics value. For small metal cylinders and 
spheres, poy =7/2In(4i/yb), po: =31b?/16, and p3= 3703 
/32. In these expressions, pg is for the electric vector 
of the incident wave parallel to the axis of the metal 
cylinder, p21 is for the electric vector perpendicular to 
the axis, 7=./—1, and Jny=0.5772. With the 6 used 
here corresponding to the critical density radius, these | 
expressions could be used to find the reflection from the 
electron clouds when X,>1 and 0<2. However, the 
values given in Table III for the maximum reflections 
from the high-density clouds are still applicable, except 
for the defocussing factor. That is, with X,>1, the 
maximum critical radius of the clouds corresponds to a 
value of b comparable to or greater than two. 

Summarizing this section, it appears that there are 
no serious errors in the approximate method for n=1, 
or for n=2 or n=3 if b>2. However, with b<2, X, <1, 
and  =2 or nm =3, polarization forces may cause wide de- 
partures from the approximate results if the radial elec- 
tron-density gradient is large. For diffuse clouds such © 
as the Gaussian models, however, such polarization 
forces are of little importance. With <2, X,>1, and 
n=2 or n=3, the approximate results are inaccurate be- 
cause the geometrical optics approximation of reflection 
from metal cylinders and spheres is no longer valid. This 
last difficulty can easily be corrected, however, by using 
the appropriate expressions for scattering from small 
metal spheres and cylinders. 


POssIBLE APPLICATIONS 


Several possible applications of the theoretical results 
will be considered very briefly. It should be emphasized 
that the method of analysis employed above is intended 
to give approximate answers for the maximum number 
of situations with the minimum amount of work. Thus, 
while some detail is lost, it is hoped that this broad view 
of radio reflections from ionized regions of fairly general 
characteristics will prove useful. 


enn A. Manning, “The strength of meteoric echoes from dense 
columns,” Jour. At. Terr. Phys., vol. 4, pp. 219-229; December, 1953. 
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The results for 7 =1 may be used to make estimates 
of layer thickness and ionization content from the 
measured behavior of radio reflections from sporadic-E 
clouds. If an auroral arc is essentially a sheet of ioniza- 
tion parallel to the geomagnetic field lines, a similar 
application of the theoretical results can be made. 

The theory as presented here for »=2 has been ap- 
plied elsewhere to describe the characteristics of radio 
reflections from meteor ionization trails.!~* The use of 
the Gaussian model has produced simple answers for 
the amplitude, duration, and polarization of meteor 
echoes which check fairly well with experiment. It has 
also been suggested that radio reflections from auroral 
rays can be analyzed as reflections from many columns 
of ionization.'° 

The results which have been obtained for »=3 may 
be applicable in explaining the complex meteor echoes 
which are observed from very dense meteor trails. For 
example, a strong reflection is obtained from a smooth 
column of ionization only when the incident and re- 
flected rays make supplementary angles with the axis 
of the column. However, the range-time displays pre- 
sented by McKinley and Millman" show a complex echo 


10 R. K. Moore, “Theory of radio scattering from the aurora,” 
Trans. I.R.E., vol. 3, pp. 217-230; August, 1952. 

1 —D. W. R. McKinley and P. M. Millman, “A phenomenological 
theory of radar echoes from meteors,” Proc. I.R.E., vol. 37, pp. 364— 
_ 375; April 1949. 
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structure with long-duration returns from other than the 
specular reflection point. The general features of these 
complex echoes would be obtained from an ionized region 
formed from point sources of ionization superimposed on 
a line source. That is, as McKinley and Millman have 
assumed, the larger meteoric particles, as they speed 
through the #-region, may produce ionization at a unt- 
form or slowly varying rate except for several localized 
regions where large ionization concentrations are cre- 
ated. This line-plus-point source model for the high- 
density meteor trails would also explain why the wave- 
length dependence for the duration of the long-enduring 
echoes, as measured by McKinley,” is nearer \4/* than 
d (see Table V). 

If the irregularities of ionization density in the 
ionosphere are regarded as being a distribution of 
ionized clouds of various sizes, shapes, and densities, 
the treatment given above for n=1, nm =2, and »=3 can 
be employed to determine the characteristics of the 
fields scattered by such ionospheric irregularities. While 
a complete description of the heterogeneous ionosphere 
is probably impossible, clues to its fine structure may 
be obtained by comparing measured results with theo- 
retical predictions based on various models. 


2D, W. R. McKinley, “Meteor echo duration and radio wave- 
length,” Can. Jour. Phys., vol. 31, pp. 1121-1135; November, 1953 
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Discussion on Optimum Patterns for Endfire Arrays” 
R. L. PRITCHARDt 


Summary—The method of synthesizing an equal-minor-lobe di- 
rectivity pattern suggested by Riblet for the broadside array was ap- 
plied recently by DuHamel to the case of the endfire array. In the 
present paper, an alternative synthesis procedure is described, based 
directly on Dolph’s method of synthesis for the broadside array. Ad- 
vantages of this alternative procedure include (a) applicability to ar- 
rays having even as well as odd numbers of elements, and (b) some- 
what simpler equations for calculating relative currents for the 
elements of the array. This alternative procedure is described here 
in detail, together with numerical examples for the seven-element 
array used by DuHamel and for a four-element array. A tabulation 
of equations for relative current amplitudes for arrays of 5, 7, 9, 11, 
and 13 elements is given in an Appendix. In conclusion, an alterna- 
tive method of overdesigning a supergain antenna is described. 


N A RECENT paper, DuHamel! has described a 
if method of synthesizing an equal-minor-lobe, or 

Tchebycheff, directivity pattern for an endfire 
linear array. This method is based on a generalization 
of the synthesis procedure suggested by Riblet,? in ex- 
tending the work of Dolph,* for the broadside array to 
include the case of arrays having an element spacing 
less than a half-wavelength. Dolph’s original method 
yields the optimum pattern only for element spacings 
greater than a half-wavelength, whereas Riblet’s meth- 
od permits an optimum pattern to be specified for any 
element spacing. However, Riblet’s method is applica- 
ble, as such, to arrays having an odd number of elements 
only, whereas Dolph’s method may be used for even or 
odd numbers of elements. 

The purpose of this communication is to point out 
that the method originally described by Dolph may be 
applied directly to the case of the endfire array. An op- 
timum pattern may be obtained for the endfire array 
for any element spacing, in contrast with the case of the 
broadside array. (It should be noted, however, that ele- 
ment spacings greater than a half-wavelength, normally 
are not used in an endfire array in order to avoid extra 
major lobes in the pattern.*) There are two advantages 
of using Dolph’s method. First, the same method may be 
applied directly to arrays having even or odd numbers 
of elements. Second, the resulting equations for deter- 
mining the relative currents for the elements (excitation 


* Original manuscript received by the PGAP, July 9, 1954. 

t 2153-A Daisy Lane, Schenectady, N. Y. 

1R. H. DuHamel, “Optimum patterns for endfire arrays,” 
Proc. I.R.E., vol. 41, pp. 652-659; May, 1953. 

2H. J. Riblet, Discussion on “A current distribution for broad- 
side arrays, which optimizes the relationship between beam width 
and side-lobe level,” Proc. I.R.E., vol. 35, pp. 489-492; May, 1947. 

§C. L. Dolph, “A current distribution for broadside arrays which 
optimizes the relationship between beam width and side-lobe level,” 
Proc. I.R.E., vol. 34, pp. 335-348; June, 1946. 

4 This was noted by DuHamel, Joc. cit., p. 655. 


coefficients) are of somewhat simpler form than the cor- 
responding equations derived by DuHamel from Riblet’s 
method. 

Application of Dolph’s method to the endfire array is 
illustrated in Fig. 1 for the five-element array used by 
DuHamel. At the right-hand side of Fig. 1 is shown the 
space factor S;(W) of the array expressed as a function 
of the auxiliary variable 


y = Bd cos ¢ — a. (1) 


Using DuHamel’s notation, ¢ is the direction angle 
measured from the axis of the array, B=27/) is the 
free-space phase constant (A=free-space wavelength), 
d is the element spacing, and a is the uniform progressive 
phase shift between adjacent elements of the array. This 
function S;(W) is essentially the same as DuHamel’s 
T2(W) shown in his Fig. 5 (except for a shift of a in the 
abscissa scale). 


z=A COS(#/2) 


Fig. 1—Graphical construction of optimum directivity pattern from 
T4(z) polynomial for five-element endfire array with d/\=3. 


In DuHamel’s method, S;(W) is determined from the 
second degree Tchebycheff polynomial 72(s) by means 
of a linear transformation of the type z=a cos w-+0. In 
the method to be described here, S;(W) is determined 
from the fourth-degree Tchebycheff polynomial T4(z) by 
means of a linear transformation of the type 


z= A cos (W/2), (2) 


as shown graphically for the example of Fig. 1. 
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For the general case of an N-element array (N even 
or odd), Sv(W) may be determined from the Tchebycheft 
polynomial 7'y_1(z), using the transformation (2). The 
constants A and a are determined by imposing the two 
conditions: 


(ajpate = 0, 2 =-20,.n-4 


I 
a 


so that | Sy(y) | 


3 
(owato— 7. so —'1 so that | Sx (y) | 


I 


where R is the specified major-lobe amplitude relative to 
unity minor-lobe amplitude, and zo y_; is related to R 
through the expression 

Ty-(Z0,n-1) = R. (4) 
Explicitly, 


Zo,v-1 = cosh {[1/(N — 1)] cosh R}, (5) 


Note that a subscript N—1 has been added to the 
parameter 2 in order to distinguish between this quan- 
tity and the parameter used in DuHamel’s method. 

By imposing (3) on (2), it follows that 


Z0,v-1 = A cos [(8d — a)/2] 
— 1=A cos [(8d + a)/2]. 


Solving these equations for A and a yields 
A = + (1/sin Bd) [(2?0,w—a = 1) + 220/n-4 Cos Bd |2/2, (6) 
a = 2tan“ {cot (84/2) [(zo,v-1 + 1)/(go,v-1 — 1)]}. (7) 


A is taken to be positive, and a normally is in the first 
or second quadrant. 

When A and a are determined, the properties of the 
directivity pattern, e.g., the width of the major lobe, 
the location of the nulls, etc., may be determined from 
_ the known properties of the Tchebycheff polynomial 
Ty-1(z) and the transformation (2). Values of the am- 
plitudes J, of the currents for the elements in the array 
may be determined from the relation 


Ty(A cos y/2) = Sy(y), (8) 


by expanding both sides of the equation in powers of 
cos (W/2), as shown originally by Dolph. In fact, the 
amplitudes of the currents for the endfire array may be 
obtained directly from Dolph’s results for J, by replac- 
ing the parameter 2 in his equations by A as defined 
above. Explicit formulas for J, for 12-, 16-, 20-, and 
24-element arrays were given by Dolph,’ but corre- 
sponding formulas for odd numbers of elements do not 
appear to have been published. Accordingly, the rela- 
tions for I, for N=5, 7, 9, 11, and 13 elements are in- 
cluded here in an Appendix. 


5 See also, D. Barbiere, “A method for calculating the current 
distribution of Tchebycheff arrays,” Proc. I.R.E., vol. 40, pp. 78-82, 
January, 1952; p. 991, August, 1952. Also, Gay. van der Maas, 
“A simplified calculation for Dolph-Tchebycheff arrays,” Jour. Appl. 
Phys., vol. 25, pp. 121-124: January, 1954. 
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EXAMPLES 


The method described above will be illustrated for 
the same example used by DuHamel of a 7-element ar- 
ray with d/\=} and with a side-lobe level of —20 db. 
From (5) above, 


Zoe = cosh [(1/6) cosh-! 10] = 1.1271. 
Substituting this result in (6) and (7) with Bd =7/2 yields 


Ai ==522 704. a = 173.15 degrees. 


This value of a is essentially the same as that calculated 
by DuHamel® except for the factor 7, i.e., DuHamel’s 
—6.8 degrees=(a—7) above. Values of the current 
magnitudes J, calculated from the appropriate equa- 
tions in the Appendix for the above value of A also are 
identical with those calculated by DuHamel.® 

An array of an even number of elements may be 
treated in exactly the same manner. For example, con- 
sider a 4-element array with d/A=} and with a side- 
lobe level of —20 db. From (5) 


Zo3 = cosh |(1/3) cosh-! 10] = 1.5408. 


(Note that this is the same value of 2) used by DuHamel 
for his 7-element array.) Substituting this result in (6) 


and (7) yields 
A = 1.8369, a = 156.0 degrees. 


The current amplitudes J, may be calculated easily 

from the explicit form of (8): 
I, = (3/2) A(A? — 1); In = (1/2)A*. 

Thus, the currents 5, =J,,¢—7@"—) ¢/2) are: 


3, = 5_1* = 6.542/78 degrees; 
Je = T_o* = 3.099/234 degrees. 


A curve of the directivity pattern (normalized space fac- 
tor) for this example is shown in Fig. 2, next page. 


Bi-DIRECTIONAL ARRAY 


The bi-directional endfire array considered by Du- 
Hamel also may be treated by the adaptation of Dolph’s 
method described here.” Eqs. (2) and (8) are applicable 
as before, but a is set =7, and A is determined from the 
condition that at ¢=0, z=20,n-1. Hence, (2) may be 
solved for 


A = 29,n-1/cos [(8d — 1)/2] = 20,n-1/sin (wd/d). (9) 


6 DuHamel, loc. cit., p. 655. : 

7 This case was described by the writer in his Ph.D. Thesis, 
“Directivity of Acoustic Linear Point Arrays,” pp. 3-48 to 3-53, 
Harvard University, October, 1950, and in Tech. Memo. No. 21, 
pp. 38-42, issued by the Acoustics Res. Lab., Harvard University, 
January 15, 1951. The material relative to equal-minor-lobe patterns 
also appears in “Optimum directivity patterns for linear point ar- 
rays,” Jour. Acous. Soc. Amer., vol. 25, pp. 890-891; September, 1953. 


42 T-R-E TRANSACTIONS—ANTENNAS AND PROPAGATION 


INDICATION OF WHEN Do.Lpu’s METHOD IS 
APPLICABLE 


For a pattern synthesis problem in general, an indica- 
tion of whether or not Dolph’s method is directly ap- 
plicable can be obtained from inspection of the re- 
quired curve of the space factor S expressed asa function 
of the auxiliary variable y=8d cos ¢ — a. If the re- 
quired curve Sy(W) for NV elements contains principal 
maxima of only one amplitude, the required Sy(y) curve 
can be obtained from the Tchebycheff polynomial 
Tw_-1(z) by a simple linear transformation of the type 
z2=A cos (W/2). This is the case for the broadside array 
when d/A=1/2, as described by Dolph? (cf. DuHamel’s 
Fig. 2) and for the endfire array as described above (cf. 
S;(W) in Fig. 1 above, or DuHamel’s Fig. 5). On the 
other hand, if the required Sy(W) curve contains princi- 
pal maxima of two different amplitudes, then it is 
necessary to use the T(z) polynomial for an array with 
an odd number 2M-+1 elements. This is the case for a 
broadside array with d/\<1/2, as shown by Riblet? 
(cf. DuHamel’s Fig. 4). The same type of curve also 
may be encountered in the case of a linear array having 
a pattern that is steered® away from the broadside di- 
rection, when d/A <1/2. 
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Fig. 2—Calculated directivity pattern for four-element endfire 
array with equal-minor-lobe excitation for d/A =}. 


ALTERNATIVE METHOD OF OvER-DESIGN FOR 
SUPERGAIN ARRAYS 


In conclusion, the writer should like to point out an 
alternative method of “over-designing” the endfire ar- 
ray in order to avoid large current amplitudes that can- 
not be adjusted with sufficient accuracy in a practical 
array to produce the desired pattern. The method de- 
scribed below may be applied also to broadside arrays. 

DuHamel suggested calculating current amplitudes 
for a larger major-lobe amplitude than is required (e.g., 
for a —20 db side-lobe requirement, design for —50 
db). This yields smaller current amplitudes at the ex- 
pense of a wider major lobe. 


_' Synthesis of a steered pattern has been discussed briefly by the 
writer (reference 7). 
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Alternatively, the pattern may be designed for the 
required minor-lobe amplitude without forcing the pat- 
tern to have as many minor lobes as is possible. For 
example, for the five-element array illustrated in Fig. 1, 
the optimum pattern has 4 minor lobes. If only 3 minor 
lobes are permitted within the visible range of y, for 
the same major-lobe amplitude, the amplitude of the 
principal maxima of S(W) (at Y=0, —2z, ...) will be re- 
duced. This will yield a smaller magnitude of currents 
for the elements, at the expense of a wider major lobe. 

The design of such a nonoptimum pattern may be 
carried out in exactly the same manner as that described 
above for the optimum pattern, except for the condi- 
tion (b) of (3). This must be replaced by the condition 
that at 6=7, z=24= —1, where zg may be selected by 
inspection of the curve of Ty_i(z) to include as many 
secondary maxima as desired, e.g., in Fig. 1, a choice of 
za= —0.7 would yield a pattern having 3 minor lobes. 
In general, the larger the value of zz chosen, the smaller 
will be the current amplitudes and the larger will be 
the beamwidth. 

For a design of this type, A and a may be determined 
as: 


A = (1/sin Bd) [(2%o,v-1 + 242) — 2%0,n-12a cos Bd]!/?, (10) 
ay == Diya { cot (Bd/2) [(Z0,—1 — 24a) /(%0,n-1 + 2a) | \ tt) 


However, za cannot be chosen completely arbitrarily. 
For each optimum pattern there is an upper limit to 2g 
(which will be defined as 24’), i.e., there is a minimum 
number of minor lobes that a pattern may have when 
“over-designed” in the manner described. In particular, 
when condition 3(a) is imposed [along with the alter- 
native 3(b)]), for ¢=0, W must not exceed 0, i.e., 
(8d —a) $0 is required. 
For the limiting value of gg=2,’, it follows that 


a= Bd, A = 20 (12) 


and that 


Za’ = 20,n-1 COS Bd. (13) 


Furthermore, the space factor S(W) in this case is ex- 
actly the same as that which may be obtained from 
Dolph’s original method of synthesis for the broadside 
array, simply by adding a uniform progressive phase 
shift of Bd radians between elements. 

In the extreme limiting case of z4’=—1, application 
of Dolph’s original method with the addition of the 
phase shift 6d actually yields the optimum pattern for 
the endfire array, which is the same as that obtained by 
the method described above or by the method of Du- 
Hamel. This case has been discussed recently by 
Rhodes.® Over-design in the manner suggested here is 
not possible for this special case. However, in general, 
it is not necessary, since the current amplitudes of such 
an optimum pattern are not large, i.e., such a pattern is 


not superdirective. 


*D. R. Rhodes, “The optimum linear array for a single main 
beam,” Proc. I.R.E., vol. 41, pp. 793-794; June, 1953. 
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To illustrate the “over-design” method described 
above, consider the example of the 7-element array, 
with d/A=%, R=10, used by DuHamel. The limiting 
value of za from (13) is z4’=0. For this value of zg, in- 
spection of the T(z) polynomial indicates that 3 minor 
lobes would occur in the pattern, instead of the 6 in the 
optimum pattern. Suppose that zg= —0.5 is chosen; this 
will yield 4 minor lobes in the pattern. Eqs. (10) and 
(11), with zo =205=1.1271, yield 


A* = 1.5204, a = 137.9 degrees. 


The current amplitudes J, calculated from equations in 
the Appendix for NV =7, with this value of A, are as fol- 
lows 


I, = 3.609, Is = 1.757, 


Note that the current amplitudes for the center ele- 
ments are approximately 15 per cent as large as those 
corresponding to the optimum pattern. This is the order 
of magnitude of reduction achieved by DuHamel for his 
over-designed optimum patterns nos. 3 and 4. Calcula- 
tion of the beamwidth of the 4-minor-lobe pattern 
yields a value of 57 degrees between the half-field in- 
tensity (—6db) points on the pattern. This, of course, 
is larger than that of the optimum pattern, but com- 
pares favorably with the calculated beamwidth for Du- 
Hamel’s over-designed optimum patterns no. 2 and 
mo. 3. 

If the limiting value of z7=0 is chosen, yielding a 3- 
minor-lobe pattern, the current amplitudes are reduced 
further (e.g., Ip=1.886), but the beamwidth is in- 
creased (e.g., 105 degrees between half-intensity points). 


APPENDIX 


CURRENT AMPLITUDES FOR LINEAR ARRAYS HAVING 
Opp NUMBERS OF ELEMENTS, WITH EQUAL- 
Mrnor-LoBE EXCITATION 


Obtained from relation: 


(N—1)/2 


Ty(A cos y/2) = Di repli COs gy: 


n=0 
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N=5 
Iy=3A'—44241 
I,;=2A4—2A? 
eh) As 
N=7 
Ip=10A'—18A*+9A?—1 
I, =(1/2)(15A*— 2444+ 9A?) 
I,=3A'—3A?4 
i,=(1/2)A* 
N=9 
Ip =35A*— 80A'+60A4—16A?-+1 
I,=28A'—60A'+40A!—8A? 
I,=14A4®—24A*+10A4 
I3=4A8—4A® 
I,=(1/2)A8 
N=11 
Ip=126A—350A*-+-350A°— 15044 25A72—1 
Te (12) [210A 1°— 5604*+525.45— 200444 254?| 
I,=60A?°— 140A8+105A®—25A4 
I3= (1/2) [45.41°—80A8+35A8| 
14=5A"—5SA8 
Ip=(1/2)A 
N=13 


Typ = 462A — 151244 1890A8— 1120A® 
+3154!—36A?+1 

1, =396A!2—1260A!0+ 15124 8— 840A ® 
+21044—18A? 

Ty= (1/2) [495A #2— 1440A 94 151248 
—672A®+105A 4} 

I3=110A!2—270A0+216A8—56A® 

In=33A¥—60A1+27A8 

I5=6A"—6A1 
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